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Abstract. We review recent progress in the construction of modified gravity models 
as alternatives to dark energy as well as the development of cosmological tests of 
gravity. Einstein’s theory of General Relativity (GR) has been tested accurately within 
the local universe i.e. the Solar System, but this leaves the possibility open that it is 
not a good description of gravity at the largest scales in the Universe. This being said, 
the standard model of cosmology assumes GR on all scales. In 1998, astronomers made 
the surprising discovery that the expansion of the Universe is accelerating, not slowing 
down. This late-time acceleration of the Universe has become the most challenging 
problem in theoretical physics. Within the framework of GR, the acceleration would 
originate from an unknown dark energy. Alternatively, it could be that there is no 
dark energy and GR itself is in error on cosmological scales. 

In this review, we first give an overview of recent developments in modified 
gravity theories including f{R) gravity, braneworld gravity, Horndeski theory and 
massive/bigravity theory. We then focus on common properties these models share, 
such as screening mechanisms they use to evade the stringent Solar System tests. Once 
armed with a theoretical knowledge of modified gravity models, we move on to discuss 
how we can test modihcations of gravity on cosmological scales. We present tests 
of gravity using linear cosmological perturbations and review the latest constraints on 
deviations from the standard AGDM model. Since screening mechanisms leave distinct 
signatures in the non-linear structure formation, we also review novel astrophysical 
tests of gravity using clusters, dwarf galaxies and stars. 

The last decade has seen a number of new constraints placed on gravity from astro- 
physical to cosmological scales. Thanks to on-going and future surveys, cosmological 
tests of gravity will enjoy another, possibly even more, exciting ten years. 
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Einstein’s theory of General Relativity (GR) has proven successful over many years of 
experimental tests [1]. These tests range from millimetre scale tests in the laboratory to 
Solar System tests and consistency with gravitational wave emission by binary pulsars. 

The standard model of cosmology assumes GR as the theory to describe gravity 
on all scales. In 1998, astronomers made the surprising discovery that the expansion 
of the Universe is accelerating, not slowing down [2, 3]. This discovery was the subject 
of the 2011 Nobel Prize in physics. This late-time acceleration of the Universe has 
become the most challenging problem in theoretical physics. Within the framework 
of GR, the acceleration would originate from an unknown dark energy. The simplest 
option is the cosmological constant, hrst introduced by Einstein. However, in order 
to explain the current acceleration of the Universe, the required value of this constant 
must be incredibly small. Particle physics predicts the existence of vacuum energy which 
provides a value for the cosmological constant, but this is typically more than 50 orders 
of magnitude larger than the observed values that assume GR. Alternatively, there could 
be no dark energy if GR itself is in error on cosmological scales. 

The standard model of cosmology is based on a huge extrapolation of our 
limited knowledge of gravity. GR has not been tested independently on galactic and 
cosmological scales. This discovery of the late time acceleration of the Universe may 
require us to revise the theory of gravity on cosmological scales and the standard model 
of cosmology based on GR. 

It is extremely timely to tackle this challenge now. Over the next hve years, 
a number of vast astronomical surveys of the galaxy distribution will be underway, 
such as the Dark Energy Survey (DES, 2012-2017) [4], the extended Baryon Oscillation 
Spectroscopic Survey (eBOSS, 2014-2018) [5], and the Mapping Nearby Galaxies at APO 
(MaNGA, 2014-2018) [6]. These surveys will dramatically transform our measurements 
of the cosmic expansion and the large scale structure of the Universe thus providing a 
new opportunity to test gravity on astrophysical and cosmological scales. Future surveys 
such as the Euclid mission (2020- )[7], the Dark Energy Spectroscopic Instrument (DESI, 
2018-) [8] and the Subaru Measurements of Images and Redshifts (SuMIRe, 2018-) [9] 
will provide an opportunity to perform ultimate tests of gravity on the largest scales in 
our Universe. 

Recently, there has been signihcant progress in developing modihed gravity theories 
that act as an alternative to dark energy (see a review [10]). Modihed gravity models 
provide interesting theoretical ideas to tackle the cosmological constant problem and 
explain the late time acceleration of the Universe. One of the challenges for modihed 
gravity models is to satisfy the stringent Solar System constraints whilst modifying 
gravity signihcantly on cosmological scales. Screening mechanisms have been developed 
to hide modihcations of gravity on small scales (see a review [11]). This research 
has been further developed into tests of GR itself via cosmological observations (see 
reviews [12, 13]). It is in principle possible to construct non-parametric consistency 
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tests of GR on cosmological scales by combining various probes of large-scale structure, 
because the Einstein equations enforce a particular relation between observables. Non¬ 
linear structure formation is complicated by screening mechanisms and new techniques 
are required to analyse it. Analytic methods and N-body simulation techniques have 
been developed to study non-linear clustering of dark matter under their influence. 
Furthermore, screening mechanisms have also opened up new possibilities to test gravity 
using astrophysical objects. 

We will review recent developments in models of modified gravity and explain 
why the realisation of a screening mechanism is important in these models. We then 
discuss how we test gravity on cosmological scales and emphasise the importance 
of understanding the non-linear structure formation. Based on the latest N-body 
simulations, we study the effects of screening mechanisms on the structures in our 
Universe. Finally we summarise the novel approach to test gravity using astrophysical 
objects. 

This review is not meant to cover all theoretical models nor observational tests 
of modihed gravity models systematically. Rather, we focus on key ideas behind 
cosmological tests of gravity whilst emphasising an interplay between theoretical physics, 
numerical astrophysics, and cosmological and astrophysical observations. For more 
thorough reviews, we recommend the reader to refer to excellent recent reviews on 
the subject Refs. [10, 11, 12, 13]. 

This review is organised as follows. In section 2, we explain the motivation to 
consider modifications of GR and why it is difficult to construct such models. Then 
we introduce several examples of models of this type. In section 3, we explain why a 
screening mechanism is important for the recently developed modihed gravity models 
and explain how it works to hide the modihcations on small scales. In section 4, we 
discuss cosmological tests of gravity. We explain how one can construct a consistency test 
of the GR model with a cosmological constant, the A Gold Dark Matter (AGDM) model 
and why it is important to develop model independent parametrisations of deviations 
from AGDM for linear cosmological perturbations. We then explain the complexities 
that arise from non-linear structure formation. We introduce analytic methods based 
on perturbation theory to describe quasi non-linear perturbations and emphasise the 
importance of non-linear physics in extracting information about modihcations of gravity 
on linear scales. We then discuss N-body simulations to study fully non-linear structures. 
In section 5, we discuss how screening mechanisms operate in dark matter simulations 
by using two representative screening mechanisms. Section 6 is devoted to novel 
astrophysical tests of gravity using clusters, galaxies and stars. We discuss the future 
outlook in section 7. 
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2.1. Problems with the standard model of cosmology 


The standard cosmology is highly successful at explaining a number of observations in 
a simple framework. It is based on two main assumptions: 


• Our universe is homogeneous and isotropic. 

• Gravity is described by General Relativity (GR). 

We use the first assumption to write down the metric for the background universe as 
the Friedman-Robertson-Walker (FRW) metric 


= -dC + a{tf ‘(^-^—+ r^dO^ j, 

and from the Einstein equations we obtain the Friedman equation 
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the Friedman equation today can be rewritten as 


1 — flrn + + Ga. (4) 

Here Hq is the present-day Hubble parameter. Many independent datasets intersect 
on the Qrn and Ga plane (see Fig. 1). However, ordinary matter occupies only 5% of 
the energy density of the Universe. We need to assume that 27% of the energy density 
of the Universe is made of dark matter and 68% of the energy density is made of the 
cosmological constant. This surprising result comes from the fact that the expansion 
of the Universe today is accelerating. The accelerating expansion of the Universe was 
found by measuring the distance to supernovae (SNe). Due to the accelerated expansion 
of the Universe, supernovae look dimmer than what we would expect in the Universe 
without a positive cosmological constant. See Ref. [14, 15] for reviews. 

The energy density associated with the cosmological constant is extremely small. 
Using the fact that the present day Hubble parameter Hq is given by Hq = 2.13h x 10“^^ 
GeV, where Hq = 100/i(km/s)/Mpc, the energy density associated with the cosmological 
constant is incredibly small 

= 8^ ~ (5) 

At the moment we do not have any compelling theory that allows such a small 
cosmological constant. The cosmological constant problem is arguably the most 
challenging problems in theoretical physics [20]. Quantum field theory predicts the 
existence of vacuum energy. Without gravity, this is the zero-point energy and it 
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Supernova Cosmology Project 
Suzuki, et al.. Ap.J. (2011) 



Figure 1. Left: Constraints on and ftm from Supernovae (SNe), Cosmic Microwave 
Background (CMB) and Baryon Acoustic Oscillations (BAO) [16]. From [16, 17] 
(published on 10 February 2012 ©AAS. Reproduced with permission). Right: The 
relative amounts of the different constituents of the Universe [18]. From [19] (©ESA 
and the Planck Collaboration). 


does not affect the dynamics. However, in GR, the vacuum energy gravitates and 
therefore the value of this vacuum energy has an important meaning. The vacuum 
energy is determined by the cut-off scale of the theory and this is typically many orders 
of magnitude larger than the cosmological constant that we need to explain the late 
time acceleration of the Universe (5). For example, for a massive particle with mass 
m, we expect that the contribution to the vacuum energy is 0{m^). Even the electron 
Tile = 0.5 MeV gives a huge contribution to the vacuum energy. It is possible to add 
a classical contribution to the cosmological constant in GR. By tuning this piece, it 
is in principle possible to realise the small cosmological constant that is required for 
the late time acceleration of the Universe. However, the problem is that this tuning is 
unstable under quantum corrections. Any new additions to the matter sector or higher 
order loop contributions spoil the tuning (see a recent review [21]). This cosmological 
constant problem was already a serious problem before the discovery of the accelerated 
expansion of the Universe, thus it is sometimes called the old cosmological constant 
problem. The discovery of the accelerated expansion of the Universe has worsened the 
problem. Now not only do we need to explain why we do not see a large cosmological 
constant but also we require an explanation for why we observe such a tiny cosmological 
constant. 

Even if we could solve the old cosmological constant problem, we still need to 
explain the origin of the accelerated expansion of the Universe. There are mainly 
two approaches to this problem. One possibility is to abandon the assumption of 
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homogeneity. Although the isotropy of the universe has been proven to high accuracy by 
the Cosmic Microwave Background (CMB) anisotropies whose temperature is isotropic 
to roughly one part in 100, 000, the homogeneity of the Universe has not been tested well. 
If the accelerated expansion of the universe is connected to the inhomogeneity of the 
Universe, this provides a compelling explanation for the problem known as “why now”. 
If the accelerated expansion of the Universe is driven by the cosmological constant, it is 
indeed quite a remarkable coincidence that flm and IIa are roughly the same order today. 
The ratio between Qm and Ha is nearly zero until recently and it approaches one rapidly 
in the future. If the acceleration is driven by the inhomogeneity caused of the growth of 
structure in our Universe, it can explain this coincidence between our existence and the 
accelerated expansion of the Universe. Unfortunately, recent studies indicate that it is 
not possible to explain the accelerated expansion of the Universe from the back-reaction 
of the structure formation in the Universe. Although, it is still possible to explain the 
acceleration by large inhomogeneities if they are not related to the structure growth. For 
example if we are living in a large void, it is in principle possible to explain the observed 
accelerated expansion of the Universe measured by SNe. However, the consensus is that 
it is difficult to explain all the observations consistently in these void models. Also 
this model requires us to abandon the Copernican principle which states that we are 
not living in a special place in the Universe. It is still an open question whether more 
sophisticated inhomogeneous models can explain the late time acceleration or not. (See 
a review [22] and references therein). 

The most popular approach is to introduce dark energy to explain the acceleration 
within the framework of GR (see a review [23]). We can replace the cosmological 
constant by dark energy with the energy density pde and the pressure Pde- In order to 
explain the accelerated expansion of the universe, the equation of state Wde = Pde/Pde 
of dark energy needs to be smaller than —1/3. The cosmological constant is a special 
case where Wde = —1- The most studied candidate for dark energy is a quintessence 
model in which the scalar held plays the role of dark energy whose energy density and 
pressure are given by 

Pde = lfi^ + V{fi), Pde=^0'-U(0), (6) 

where 0 is a scalar held. If the scalar held’s kinetic energy is sub-dominant compared 
to its potential energy, we can realise the negative pressure that is required to explain 
the accelerated expansion of the Universe. It is also possible to hnd a tracking solution 
where the scalar held energy density follows that of matter Pde/Pm =const. for an 
exponential potential V oc exp(—A(/)). This means that the hnal density of the dark 
energy is insensitive to the initial conditions, which could solve the why now problem. 
However, for this simple example, the scalar held never dominates the Universe and 
it cannot explain the accelerated expansion of the Universe. In order to achieve the 
acceleration, we need to add another exponential potential and the time when dark 
energy dominates over the matter energy density is controlled by the parameters in the 
potential. This is exactly the same situation in the ACDM model and we can therefore 
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no longer solve the why now problem. Another difficulty of the quintessence models 
is that the parameters in the potential are unstable under radiative corrections. In 
order for the scalar held to be a successful candidate for dark energy, the mass of the 
scalar held needs to be of the order Hq and this mass is not protected under quantum 
corrections. Due to these theoretical difficulties, dark energy is often treated as a huid 
with a given equation of state as a function of scale factor tCde(a) without specifying its 
physical origin. 

2.2. Why modified gravity? 

The cosmological constant problem is of fundamental origin thus it is important to 
reconsider all the assumptions that we make in the standard model of cosmology. One of 
the assumptions is that gravity is described by GR on all scales. However we should bear 
in mind that we tested gravity only in our local universe. There is an interesting example 
in history of where unexpected observations lead to a revolution in gravitational physics. 
When the precession of perihelion of Mercury was found, a French mathematician, 
Urbain Le Verrier hypothesised that there existed a planet called Vulcan between the 
Sun and Mercury that caused the anomaly. In fact he succeeded in predicting the 
existence of Neptune using the same technique. Vulcan was never found and we now 
know that it was Newton’s gravity that needed to be revised by Einstein. This story 
tells us that sometimes we need to be open-minded when we encounter unexpected 
observations. 

Fig. 2 summarises where we have tested gravity [24]. Gravity is parametrised by 
the two quantities, the gravitation potential GM/r and the curvature of space GM/r^ 
for a spherical object with mass M and the radius r. GR is very well tested in the 
Solar System and also by binary pulsars. However it is especially not well tested in the 
low curvature regime. We need dark matter to explain the rotation curves of galaxies, 
clusters of galaxies and the formation of large scale structure of the Universe. There 
have been many attempts to explain dark matter by modihcations of gravity (see [25] 
for a review). In this review, we do not cover these developments. Given the myriad 
of observations, it is becoming almost impossible to explain all the observed evidence 
for dark matter by modihcations of gravity. In the next decades, progress in direct and 
indirect detection experiments will shed light on the origin of dark matter (see a review 
[26]). In this review, we assume that dark matter is the dominant constituent of the 
Universe at late times. At the lowest curvature, we require dark energy to explain the 
accelerated expansion of the Universe. Here only cosmology provides a means to test 
gravity. 

Modihcations to gravity could provide an interesting solution to the cosmological 
constant problem. The hrst idea is to modify the way in which gravity responds to 
the cosmological constant [27]. An example is provided by a braneworld model in 6D 
(see [28, 29] for reviews). Braneworld theories postulate that we are living in a 4D 
membrane embedded in a higher dimensional spacetime (see [30] for a review). In 6D, 
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Figure 2. Left: A parameter space for gravitational fields. Right: The experimental 
version of the parameter space. See Ref. [24] for details. The horizontal lines in the left 
figure indicate the background curvature of the Universe at Big Bang Nucleosynthesis 
(BBN) and Last scattering, and the curvature associated with A. Some of the 
label abbreviations are: SS=planets of the Solar System MS=Main Sequence stars, 
WD=white dwarfs, PRSs=binary pulsars, NS=Neutron stars, BH= stellar mass black 
holes, MW=the Milky Way, SMBH=supermassive black holes. PPN= Parameterised 
Post-Newtonian regime, Inv.Sq.=laboratory tests of the inverse square law of the 
gravitational force, Atom=atom interferometry experiments, EHT=Event Horizon 
Telescope, ELT=the Extremely Large Telescope, DEET4=a hypothetical stage 4 
experiment of dark energy, Eacility=a futuristic large radio telescope such as the Square 
Kilometer Array. From [24] (published on 20 March 2015 ©AAS. Reproduced with 
permission). 


there is an interesting property of gravity that canses the cosmological constant on a 
brane not to gravitate bnt merely cnrve the extra two dimensions. Thns even if there is 
a hnge contribution to the cosmological constant in the 4D brane where standard model 
particles live, the 4D spacetime remains flat. In this way we can avoid the cosmological 
constant problem. Although it is difficult to hnd a concrete model that realises this idea 
in practice, this provides an interesting insight into a possibility that modihcations of 
gravity might address the fundamental problem of the cosmological constant. 

Modihed gravity models also provide interesting ideas to explain the late time 
acceleration of the Universe. One possibility to explain the acceleration is to add a 
tiny mass to the graviton (see [31] for a review). The smallness of the mass can be 
protected by the restoration of a symmetry. A mass term breaks the diffeomorphism 
invariance that is present in GR in the massless limit. Due to this restoring of the 
symmetry in the massless limit, it is technically natural to have a small mass. More 
precisely if we compute the quantum loop corrections to the mass, the mass is shifted 
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by —)■ because the quantum corrections arise only from the mass term 

that breaks the symmetry. Thus it is natural to tune the mass to be small to account 
for the late time acceleration ~ Hq where Hq is the present day Hubble scale as this 
tuning is not spoiled by the quantum corrections unlike the cosmological constant. 

Another interesting idea is self-acceleration. Even without the cosmological 
constant, the expansion of the Universe can accelerate due to the modihcation of gravity. 
A typical example is again provided by a braneworld model. In this case, a 4D brane 
is embedded in a 5D spacetime. Gravity in the 5D bulk spacetime is described by 
the 5D Einstein gravity [32]. On the 4D brane there is an induced gravity described 
by the 4D Einstein-Hilbert action. Moreover, there is a solution where the expansion 
of the Universe is determined by the ratio between the 4D and 5D Newton constant 
H oc G 4 /G 5 without the cosmological constant [33]. 

Unfortunately, we still do not have a model that realises these new ideas in a 
consistent way. However, clearly it is worthwhile challenging GR in order to address 
the fundamental problems in theoretical physics. Modihed gravity models also provided 
motivation for developing cosmological tests of gravity (see reviews [12, 13, 11]). It is the 
right time to challenge GR on cosmological scales. As mentioned in the introduction, 
there will be a number of astronomical surveys aiming to reveal the nature of dark 
energy. Although it is possible to construct model independent tests of gravity on large 
scales, it is important to study theoretical models to know what kind of deviations from 
GR one should search for. Also in order to exploit the vast information available on 
non-linear scales, we need concrete theoretical models. 

2.3. Why is it difficult to modify GR? 

We immediately face many problems when we modify GR. There is a powerful theorem 
known as Lovelock’s theorem. Lovelock’s theorem proves that Einstein’s equations are 
the only second-order local equations of motion for a metric derivable from the action 
in 4D. This indicates that if we modify GR, we need to have one or more of these: 

• Extra degrees of freedom 

• Higher derivatives 

• Higher dimensional spacetime 

• Non-locality 

Once we introduce these additional ingredients into the theory beyond GR, we need to 
check the theoretical consistency of the model. First we need to check that the solutions 
are stable. There are several kinds of instabilities which can be illustrated by a simple 
scalar held example described by the action 

S'= y dtd^x{Ktfi^ — Kx{di(j)){d'^(j)) — (7) 

The tachyonic instability arises when the scalar held has a negative mass squared < 0. 
This instability is not necessarily catastrophic if the instability time scale \m\~^ is long 
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enough. Another instability arises when the gradient term has a wrong sign < 0. 
This instability is more severe than the tachyonic instability as the instability time scale 
is determined by the wavenumber of the mode that we are interested in, so the time scale 
becomes shorter and shorter on small scales. On the other hand, the ghost instability, 
which arises when the time kinetic term of the scalar held has a wrong sign Kt < 0 , is 
more severe. At the quantum level, the vacuum is unstable and decays instantaneously. 
To avoid the instability, it is required to introduce a non-Lorentz-invariant cut-off in the 
theory (see [34] for a review). 

Another common problem is known as the strong coupling problem. In addition to 
the kinetic term, the scalar held can have non-linear interaction terms. An example of 
the non-linear term that appears later is 

>S'non—linear d 31-^n0((90) . (8) 

This non-linear interaction becomes important at energy scale higher than A 3 and 
quantum corrections generate terms that are suppressed by A 3 . Then in general we 
loose control of the theory beyond A 3 . This strong coupling scale is often associated 
with the energy scale related to the accelerated expansion of the Universe, Hq, which is 
extremely small compared with the scale of gravity Mp\. Thus we will often hnd that 
the strong coupling scale is rather low in modihed gravity models. This means that we 
need to treat these theories as an effective theory, which is valid only at energy scales 
lower than A 3 . 

Once the theory satishes the requirements for theoretical consistency, it also needs 
to pass observational tests. First of all, the theory needs to satisfy the stringent Solar 
System constraints. The deflection angle 9 of stars due to the Sun is observed to be [35] 

9 = (0.99992 ± 0.00023) x 1.75", (9) 

where 1.75" is the prediction of GR. Another prediction of GR is time delation due to 
the effect of the Sun’s gravitational held. This was measured very accurately using the 
Gassini satellite as [36] 

At = (1.00001 ± 0.00001)AtGK. (10) 

Any modihed theory of gravity needs to satisfy these stringent constraints on deviations 
from GR in the Solar System. On cosmological scales, observations become better 
and better and at the moment, the standard AGDM model survives all these improved 
measurements. In the background, the expansion of the Universe should look very 
similar to that of the AGDM model. The equation of state for dark energy is constrained 
to be [37] 

Wde = -1.006 ± 0.045. 


Modihed gravity models need to pass all these observational tests. 


( 11 ) 
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In this section, we discuss several representative modified gravity models that have been 
studied recently. As we mentioned in the introduction, the purpose of this review is not 
to cover all possible modified gravity models in a systematic way. We only review models 
that are relevant for later discussions. For a more complete review of various modified 
gravity models, see Ref [10]. 


2 . 4 . 1 . Brans-Dicke gravity Let’s first start from the simplest modified gravity model, 
Brans-Dicke gravity [38]. This theory includes a scalar field non-minimally coupled to 
gravity 


S = 


d^Xy/^(^fiR — + j d^X^/^Cr, 


( 12 ) 


We consider a non-relativistic source Tg = 


p. Using the quasi-static approximation 
to ignore time derivatives of perturbations compared with spatial derivatives, the 
perturbations of the metric 

d-s^ = —(1 -I- 2^>)dt^ + (1 ~ 2^)5ijdxMxfi (13) 

and the scalar field perturbation ijj = i/jq + p obey the following equations 


= AttGp - -V>, 


(14) 

(15) 

( 16 ) 


(3 -h 2ujbd)'^‘^T = -SttGp, 

$ - T = 99. 

Note that the scalar field that is non-minimally coupled to gravity gives an effective 
anisotropic stress through its perturbations, modifying the relation between <h and T. 
These equations can be rewritten as 

= AttGpp, T = (17) 

where 

4: + 2ubd 1 + 

3 + 2ubd' 2 + ubd' ^ ^ 

We recover GR in the large ubd limit. Indeed, imposing the Solar System constraints 
Eqs. (9) and (10), we obtain |?7 — 1| = (2.1 ± 2.3) x 10“^. Then the constraint on ubd is 
given by ubd > 40, 000. Once we impose this constraint on the parameter of the model 
cobd: the theory is basically indistinguishable from GR on all scales. This is one of the 
main problems of modified gravity models. The Solar System tests are so stringent that 
if they are imposed on the parameters of the model, there is virtually no room to modify 
gravity on cosmological scales. 


2 . 4 . 2 . f{R) gravity Another popular model is f{R) gravity where the Einstein-Hilbert 
action is generalised to be a function of the Ricci curvature [39, 40] 


S = / (fxX^jf{R) + f d*Xs/^C„. 


(19) 
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See [41, 42, 43] for a review. The equation of motion for the metric is fourth order thus 
this model can be classihed as a higher derivative theory. However, it is possible to 
introduce a scalar held and make the equation of motion second order. The action is 
equivalent to 

S = f ifW + (fl - 4>)f{4>)) ■ (20) 

By taking a variation with respect to 0, we obtain {R — = 0. As long as 

f"{4>) 0, and i? = 0, we recover the original action. By dehning ip = f'{(p) and 

V = f{(p) — (pf'{(p), the action can be written the same as the one for Brans-Dicke 
gravity with a potential 

S = J d^xV^{iPR-V{iP)y ( 21 ) 

Comparing this with the action (12), we notice that the BD parameter is given by 
^BD = 0. As such, if we ignore the potential, this model is already excluded by the 
Solar System constraints. However, by choosing the potential, i.e. the form of the f{R) 
function appropriately, it is possible to incorporate a screening mechanism known as 
the chameleon mechanism to evade the Solar System constraint as we will see later 
[44, 45, 46, 47]. 

In general, the scalar tensor theory is described by the following action 

S = j d‘^x^/^{^R - {dipf -V{ip)^ + j d'^xy/^Cm{gi,u)-{‘2^‘2) 

By a conformal transformation = A{(pYg^i, and a redehnition of the scalar held, we 
can transform the action to the Einstein frame 

S = j d‘^x^/^(^R-]^{d4>f -V{(t))^ + J d^x^^Cm(^A{(pfgf,^y (23) 

In this frame, the scalar held is directly coupled to matter. 

Like the quintessence models, a choice of the function f{R) leads to many models 
but the successful models for the late time cosmology share the same features [46, 48, 49]. 
In the high curvature limit, the function looks like 

^n+l 

f(R) = R-2A+\fno\^, (24) 

where R is the curvature today. This model requires an ehective cosmological constant 
to explain the observed accelerated expansion of the universe. However, it is possible 
to hnd a function f{R) so that this constant disappears in the low curvature limit. The 
correction to the ACDM disappears in the high curvature limit i? 3> .R. As we will see 
later, the Solar System constraint imposes the condition Ifml < 10“®. The background 
cosmology is indistinguishable from ACDM if |//jo| is small. In the Einstein frame, the 
potential and the coupling function take the following form 

— A(0)2 = eV^7/Mp,^ 

A/pi } 


(25) 
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and the scalar field in the original Jordan frame if = ffi = df /dR is related to as 


4 > =-y ^Mpi\og{l + fn). ( 

Cosmology of f{R) models have been studied intensively. See references in [42, 43]. 


2.4-3. Braneworld gravity The idea of braneworld was inspired by the discovery of 
D-branes in string theory (see [30] for a review). In this model, we are living on a 
4D membrane (brane) in a higher dimensional spacetime (bulk). The standard model 
particles are confined to the brane while gravity can propagate throughout the whole 
spacetime. The simplest model known as the Dvali-Gabadadze-Porrati (DGP) model is 
a 5D model described by the action [32] 


.Ml 

2 


df'xM-T>)g‘d^lR + ^ I d^Xy/^(^^‘^^R +Cm^, 


where ^^^R and are 5D and 4D Ricci curvature respectively. The ratio between 
the 5D and 4D Newton constant Vc = M|/2M| is called the cross-over scale and it is a 
parameter of the model. This model provides an interesting example of self-acceleration. 
The Friedman equation in this model is given by [33] 

H^ = ±- + —p. ( 28 ) 

At early times Hrc 3> 1, we recover the usual 4D Friedmann equation. On the other 
hand, at late times, the Hubble parameter approaches a constant H —)■ 1/rc in the 
upper branch of the solution. Thus the expansion of the Universe accelerates without 
the cosmological constant. This is known as the self-accelerating branch. On the other 
hand, the lower branch solution, often called the normal branch solution requires the 
cosmological constant to realise the accelerated expansion of the Universe. In order to 
recover the standard cosmology at early times, the cross-over scale needs to be tuned as 
hM so that the modification of gravity appears only at late times. See references 
in [30] for cosmological studies of the model. 

If we study the perturbations around this background under the quasi-static 
approximations, we find that gravity is described by the BD theory with an additional 
non-linear interaction term [50, 51, 52] 

= AnGa^p - ^VV, (29) 

2 

{‘i + 2uBD{a))V‘^(p+'^ -{ViVjipf =-87rGa'^p, (30) 

(XL 

^ - A/ = ip, (31) 

where the BD parameter ujbd{o) is given by 

ubd = -1), fi{a) = l-2rcH[l + (32) 

The scalar field originates from the bending of the brane in the 5D bulk spacetime. 
Note that the BD parameter is always negative in the self-accelerated branch and the 
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scalar field is mediating a repulsive force. This is the manifestation of the problem that 
this solution suffers from a ghost instability [53, 54, 55, 56, 57, 58]. In addition, at late 
times ubd ~ 0(1) thus by linearising the equation, we again hnd that this theory would 
be excluded by the Solar System constraints. However, the coefficient of the non-linear 
interaction term is very large, ~ This non-linear term becomes important 

even if gravity is weak. This is responsible for the screening mechanism known as the 
Vainshtein mechanism as we will see later. Interestingly, the Solar System constraints 
impose the constraint Tc > 0.177,]"^. This constraint is independent of the requirement 
to recover standard cosmology at early times Tc ~ . Thus any improvement on this 

constraint will give severe constraints on the model as we will see in section 3.5. 

An interesting extension of the DGP model is provided by 6D models. As we 
explained in section 2.2, 6D spacetime has the property that the cosmological constant 
on a 4D brane does not gravitate. However, the singular structure of the brane is far 
more complicated compared with the 5D case. There have been many attempts to 
extend the 5D DGP to 6D [59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71]. 


2.5. Galileons, Horndeski theory and beyond 


The non-linear interaction that appears in the DGP braneworld model has a very special 
property [53, 72]. The equation of motion in the flat spacetime can be derived from the 
action (8). It looks like this action leads to an equation of motion that contains higher 
order derivatives, but this is not the case as can been seen from the fact that the equation 
of motion (30) contains only second order derivatives. Furthermore,this action has the 
following symmetry in held space 


diifi -|- c^. 


(33) 


where is a constant vector. Due to its similarity with the Galilean symmetry in 
Newtonian gravity, the scalar held described by this action is called the galileon [54]. 
Interestingly, in 4D spacetime, there are only three Lagrangians, in addition to the 
canonical kinetic term, that lead to second order equations of motion with the galileon 
symmetry; 

Cf = (34) 

£f' = [(□•(>)" - . (35) 

£f' = -j(S^)"[(n^)’ - 3n4>(8,d,.l>f + 2(d,a,.j>n (36) 

There are many other expressions for the same action that are related by integrations 
by parts. The most convenient form to understand why the equation of motion contains 
only second derivatives can be constructed by making use of the Levi — Civita epsilon 
tensor to write the Lagrangian for the galileons in a compact form. Using the following 
property: 


= -(4-n)!n!ilf‘-W, 


7l...74_nQ;i...Q:n 


.7l...74-n/5l.../3n 


(37) 
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where the square brackets represent normalised anti-symmetric permutations, we can 
write the galileon Lagrangians as: 

:= ^(4)0102(034), 

£|al ^ ^(6)0102(034056), 

£|al ^ ^(8)0102(034056078), 

(38) 


Where we have dehned = .. ^ gi 35 ...i 2 ii 22 -i 24 -ng 246 ... which has been 

2n (4—n)! uii/2---i^4-n 

written in short hand as £( 2 n) and the numbers are short hand for labeled indices: 
{/ii/i 2 • • •}• Furthermore, we have that 0^i..,;i„ = d^^... 

With this notation it is very easy to see that the variation of these Lagrangians would 
never have higher than two derivatives. For instance, taking the variation of £5 gives 
us: 


0 = 5^5 = y dfxdCf 

= y'd"^a;£(8) 250102(034056078) + 30102 (<5034056078) 

= j d"^a;£(8) — 2(9i (02034056078) ~ 30304(0012056078) | <^0 
= —5 f d^X £(8) (012034056078)- 


(39) 


Where we have integrated by parts and found that the only term to survive the 
summation with the totally antisymmetric tensor £( 8 ) has, indeed, only derivatives of 
second order. 

As shown by Ref. [54], these galileon interactions can be constructed as the short- 
distance limit of a conformally invariant scalar held Lagrangian. Due to the conformal 
invariance, we generally expect that the maximally symmetric solutions including the 
self-accelerating de Sitter solution to describe the vacuum. These conformal invariant 
theories can be written in terms of a particular combination of the curvature invariants 
built out of an effective metric [54, 73]. These galileon interaction terms 

can be also constructed from the action of a probe brane hoating in a higher dimensional 
spacetime in the limit in which only a scalar held degree of freedom remains [74]. If 
we consider a probe brane in a 5D Anti-de Sitter spacetime, the resultant galileon 
action was shown to be equivalent to the conformal invariant theories constructed from 
g^J.u = [73]. Also it is possible to construct galileon interactions for multiple scalar 

helds. See [75, 76, 77] for reviews and references therein. The galileon interactions can 
be also extended to vectors [78, 79, 80]. 

In curved space, if we covariantise the action naively by replacing the partial 
derivatives by covariant derivatives, the equation of motion now contains the third 
derivative of the metric for the quartic (£ 4 ) and quintic (£ 5 ) galileons [81]. In general, 
higher derivative theories lead to the so-called Ostragradsky ghost. f{B) gravity avoids 
this instability as the equation of motion can be rewritten in the second order form 
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by introducing a scalar field. In fact, it has been claimed that in the case of a naive 
covariantisation of Eq. (38), the equations of motion can be written in the second order 
form despite the appearance of the higher derivative terms due to a hidden constraint 
[82] (see also [83, 84]). Here we first discuss a way to remove the third derivatives 
by introducing a counter term. The covariant action that leads to the second order 


equations of motion for metric and the scalar field is given by [85] 

A = (40) 

A = ^(V0)"i? - [{Dfif - (V^V.0)2] , (41) 

o 

- 3n0(V^V.0)2 + 2(V^V.0)3]. (42) 

The non-minimal coupling terms for £4 and £5 are the counter terms that are necessary 
to remove higher derivatives in the equations of motion. 

These actions can be generalised further, leading to the Horndeski action [81]: 

C2 = K{fi,X), (43) 

Cs = -G3{fi,X)Dfi, (44) 

£4 = Gfifi,X)R + G,x{fi,X) [(□0)2 - (V^V.0)2] , (45) 

£5 = G' 5 ( 0 ,X)G^"V^V ,0 

-^G,x{fi,X)[{Dfif - 3d0(V^V.0)2 + (46) 


where X = —d^d^(j)/2 and K, G3, G4 and G5 are free function of the scalar field and X. 
This action was originally found by Horndeski in 1974 [ 86 ]. The Horndeski action was 
revisited by Ref. [87] and it was shown to be equivalent to the action (46) by Ref. [ 88 ]. 

Expanding the metric and scalar field around the Minkowski background, we obtain 
the galileon interactions. In addition, there appear additional couplings between tensors 
and scalars due the non-minimal coupling [89, 90]. We consider a Minkowski background 
with a constant scalar field 0 = 0 o and study the deviations around it, namely 

0 = 00 TT, . (47) 

Notice that for the background to be a solution, one needs K = dK/d(j) = 0 at 0 = 0o. 
We expand the Horndeski action (46) in terms of the fluctuations (47), using the 
following assumptions: the fields tt and are small, hence we neglect higher order 
interactions containing the tensor perturbations as well as terms containing higher 
order powers of the scalar fluctuation tt and its first derivatives only. On the other 
hand, we keep all terms with second order derivatives of tt that preserve the galileon 
symmetry (33). The kinetic term for metric perturbations is contained in Gfi(l),x)R. 
Thus we introduce the Planck scale as G 4 = Mpj/2 and canonically normalise as 
h^i, = Mpihfj,iy. We further introduce a new mass dimension A 3 so that terms involving 
second order derivatives of tt have the appropriate dimensions. By explicitly expanding 
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the Horndeski action (46) under these prescriptions, we hnd [90] 


C‘« = - -hr£-\Xf, 


+ Can + Ccf + Ccf' + 


Ar* 


Ar‘ 




A9~5 


where are the galileon terms (36) and the tensor-scalar couplings Xjfl are given by 


fii/i 
(0 


(48) 


vw = -win] - 

x'^J = nj, - |n]n„„ + - [n^]). 

V'? = 6nj„ - enjjn] + 3n„„(|n]^ - [n^]) 

-rw(|n]’-3in]in"] + 2|n=]), (49) 


where the shorthand notations 11 ^,^ = d^d^^TT, 11 ”^, = and [n"-] = 11 ”^^ are 

introduced and is the linearised Einstein tensor. We dehne seven dimensionless 

parameters f, rj, /i, u, zu, a, (3 as follows [89] 

= M-pfi, Kx — 2G30 = rj, —Gsx + “iG^^x = 


Ai^ 


A;/pi 2 u 

G 4 X — G^^ = G 4 XX — = jj-, 


G 5 X — 


3Mpi 3w 

-p, ^bXX — 


Ai 


(50) 


where all functions are evaluated at the background, 0 = 00 and X = 0, and 
6*30 = dG'i/dfi^G'ix = dGs/dX etc. 

The above action respects the galileon symmetry (33). The coupling between the 
metric and scalar perturbations h^’^xj^J and h^^xj^J can be eliminated by the local held 
redehnition [91] 


h^u h^i, ‘2^xrj^i/ , 


(51) 


SO that the action (48) now becomes 


£=« = + ’iPXnan 


ji + 6af ^gai I iz + 2q;^ -|- A/3^ I zu + 1Oq; 0 ^g^i 

» Cl -^3 + » R + » O '^,'1 


Ai 


2 A 6 




A^ 

2 Mpi 


A 1 

ttT + 


Ml 


PI 


A/piAi " " 


(52) 


However, the coupling h^'^X^J cannot be removed by a local held redehnition [91]. 
The transformation (51) introduces a coupling between tt and via the de-mixing of 
scalar from gravity. This provides the ehective action that describes the scalar-tensor 
interactions that satishes the galileon symmetry. 

Although we started from the Horndeski action, the non-linear completion of this 
ehective action is not necessarily the Horndeski action. As we have already seen, the 
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cubic galileon term appears in the decoupling limit of the 5D braneworld model. Massive 
gravity is another example with a decoupling limit described by this effective action, 
as we will see later. Another simple example is given by an extension of the Horndeski 
action called beyond Horndeski [82]. 

In order to introduce beyond Horndeski gravity, first we formulate the Horndeski 
theory in the unitary gauge [82, 92]. In the unitary gauge, the scalar held is only the 
function of time 0 = 0(f). We perform the 3+1 ADM decomposition of space-time 
described by the line element 

ds^ = g^ydx^dx'' = —N'^dt^ + hij{dx'' + NMt){dx^ + NMt ), (53) 

We introduce the extrinsic curvature dehned by 

= h^n^-x, (54) 

where = {—N, 0, 0, 0) is a unit vector orthogonal to the constant t hyper-surfaces 
Tt and a semicolon represents a covariant derivative. We also introduce the three- 
dimensional Ricci tensor on Then we can construct a number of 

geometric scalar quantities: 

K = K%, S = A = U = (55) 

The Lagrangian that leads to second order equations of motion can be expressed in 
terms of the geometric scalars introduced above as [82, 92] 

L = AfiN, t) + AfiN, t)K + AfiN, t){K^ - S) + BfiN, t)n 
+A,{Nfi)K^ + {U - Kn/2) , (56) 

where K^ = K^ - + 2 K^^K^^^K^x■ 

By expressing the Horndeski action using the geometric scalars in the unitary gauge, 
we can hnd the relation between the coefficients Gi in Eq. (46) and Aj, Bi in Eq. (56) 

= 6-2 - XF 3 + , As = 2{-Xf/^F:,,x - 2v^G4,0 , 

A 4 = —G 4 + 2 XG4 ^^x + AG5^(^/2 , H 4 = G 4 + XiG^^if, — F^^fi)/2 , 

As = -{-Xf/^G^^x /^, ^5 = -V^F ,, (57) 

where +3 and +5 are auxiliary functions obeying the relations G 3 = +3 + 2 XF‘i^x 
and G^ x = + F^^x- Since X = —<j)'^(t)/N‘^ in unitary gauge, the functional 

dependence of 0 and X can translate to that of t and N. 

The Horndeski theories satisfy the following two conditions 

A 4 = 2 XB 4 X — B 4 , As = —-XB^ x ■ (58) 

It is possible to go beyond the Horndeski theory without imposing the two conditions 
(58) [82]. This generally gives rise to derivatives higher than second order, but it does 
not necessarily mean that an extra propagating degree of freedom is present [93, 94]. 
Interesting, if we naively covariantise the galileon action of the form given by Eq. (38), 
this is included in the beyond Horndeski theory. 
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The Horndeski theory does not specify how matter couples to gravity. Another 
possibility for beyond Horndeski theories is to consider a metric 

X)g^^ + r(0, X)d^(j)du(j), (59) 


and couple matter to gravity using g^^, [95, 83, 84], The disformal coupling without the 
kinetic term dependence r(0) has attracted a lot of interest recently (see for example 
[96] and references therein) whereas the kinetically dependent disformal coupling has 
been used to show the connection between beyond Horndeski theories given by Eq. (56) 
and the Horndeski theory [92], 

Given the number of free functions in these theories it is difficult to discuss 
cosmology in a general way. One interesting feature of this model is that it is possible 
to explain the accelerated expansion of the Universe without introducing a potential. A 
simple example is given by the Horndeski theory with [97, 98] 

A'(A') = -A', G3(A') = (jfix'j . (60) 

Note that the sign for the kinetic term is opposite to the normal scalar held. However, 
this does not imply that the perturbations are ghostly when expanding around the non¬ 
trivial cosmological background. This theory admits the tracker solution for the scalar 
held given by 


3GxH' 


and this leads to the Friedmann equation of the form 
2 




I 


H 

Ho 


_ ^ 

' 2n-l 




.-3 


(61) 


(62) 


This Friedman equation is similar to the one that we found in the DGP braneworld 
which realised self-acceleration. The covariant galileon model described by Eq. (40) 
corresponds to n = 1. The perturbations around this solution is stable despite the 
wrong sign for the kinetic term thanks to the cubic galileon term. See Refs. [99, 100, 
101, 102, 103, 104, 105, 106, 107, 108, 109, 97, 110, 111, 112, 113, 114, 115, 116, 117, 
118, 119, 120, 121] for an incomplete list of studies of cosmology in this class of models. 
Galileons have also been used to construct general models of inhation. See references in 
[11, 77, 122]. 


2.6. Massive gravity/bigravity 

Another theory that is closely related to galileons is massive gravity (see [31] for a 
review). This section is based largely on a review Ref. [123]. Gan the graviton have 
a mass? Attempts to answer this question date back to the work by Fierz and Pauli 
(FP) in 1939 [124]. They considered a mass term for linear gravitational perturbations, 
which is uniquely determined by requiring the absence of ghost degrees of freedom. The 
mass term breaks the gauge invariance of GR, leading to a graviton with Eve degrees 
of freedom instead of the two found in GR. There have been intensive studies into 
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what happens beyond the linearised theory of FP. In 1972, Boulware and Deser fonnd 
a scalar ghost mode at the non-linear level, the so called sixth degree of freedom in 
the FP theory [125]. This issne has been re-examined using an effective held theory 
approach [126]. They introduced scalar helds known as Stiickelberg helds to restore 
gauge invariance so that the theory is manifestly convariant. In this formulation, the 
Stiickelberg helds describe the additional degrees of freedom in massive graviton. They 
acquire non-linear interactions containing more than two time derivatives, signalling the 
existence of a ghost. In order to construct a consistent theory, non-linear terms should 
be added to the FP model, which are tuned so that they remove the ghost order by 
order in perturbation theory. 

Interestingly, this approach sheds light on another famous problem with FP massive 
gravity; due to contributions from the scalar degree of freedom, solutions in the FP 
model do not continuously connect to solutions in GR, even in the limit of zero graviton 
mass. This is known as the van Dam, Veltman, and Zakharov (vDVZ) discontinuity 
[127, 128]. Observations such as light bending in the Solar System would exclude the FP 
theory, no matter how small the graviton mass is. In 1972, Vainshtein [129] proposed 
a mechanism to avoid this conclusion; in the small mass limit, the scalar degree of 
freedom becomes strongly coupled and the linearised FP theory is no longer reliable. In 
this regime, higher order interactions, which are introduced to remove the ghost degree 
of freedom, should shield the scalar interaction and recover GR on sufficiently small 
scales realising the Vainshtein mechanism. 

In order to avoid the presence of a ghost, interactions have to be chosen in such a 
way that the equations of motion for the scalar degrees of freedom contain no more than 
two time derivatives, i.e the galileon terms. Therefore, one expects that any consistent 
non-linear completion of FP contains these galileon terms in the limit in which the scalar 
mode decouples from the tensor modes, the so-called decoupling limit. This turns out to 
be a powerful criteria for building higher order interactions with the desired properties. 
Indeed, following this route, de Rham and Gabadadze constructed a family of ghost-free 
extensions to the FP theory, which reduce to the galileon terms in the decoupling limit 
[130]. Adopting an effective held theory approach, the basic building block is a tensor 
corresponding to the covariantization of metric perturbations, namely: 

9liu T hfiu — = 0^(f> (63) 

The Stiickelberg helds 0“ are introduced to restore reparameterisation invariance, hence 
transforming as scalar from the point of view of the physical metric [126]. The internal 
metric rjag corresponds to a non-dynamical reference metric, usually assumed to be 
Minkowski space-time. Therefore, around hat space, we can rewrite as 

+ ffiSyd^Tr^ + (64) 

From now on, indices are raised/lowered with the dynamical metric gfj_y, unless otherwise 
stated. For example, 77^ = g^'^Hpy. Moreover, the Lagrangian constructed from H^y 
is invariant under coordinate transformations ^ provided tt^ transforms as 

TT^ —)■ TT^ -|- (65) 
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The dynamics of the Stiickelberg helds 0" are at the origin of the two featnres: 
the BD ghost excitation and the vDVZ discontinnity. With respect to the hrst issne, 
as noticed by Fierz and Pauli, one can remove the ghost excitation, to linear order in 
perturbations, by choosing the quadratic structure — H"^. When expressed in the 

Stiickelberg held language, one of the four Stiickelberg helds becomes non-dynamical. 
However, when going beyond linear order, this constraint disappears, signalling the 
emergence of an additional ghost mode [126]. Remarkably, Ref. [130] has shown how to 
construct a potential, tuned at each order in powers of to hold the constraint and 
remove one of the Stiickelberg helds. Even though the potential is expressed in terms 
of an inhnite series of terms for it can be re-summed into the following hnite form 
[131] 

U =-rn^ \U2 + , (66) 

where q;„ are free dimensionless parameters, = n!det„(/C) and the tensor )Cfi is 
dehned as 

K; = v-(k<r^)l (67) 

The square root is formally understood as The determinant can be 

written in terms of traces as 

det(/C) = (tr ICfi — tr (/C^), 

det(/C) = (tr/C)^ — 3(tr/C)(tr/C^) -|- 2tr/C^, 

det(/C) = (tr/C)^ — 6 (tr/C)^(tr/C^) -|- 8 (tr/C)(tr/C^) -|- 3(tr/C^)^ — 6 tr/C^ . 

All terms det„(/C) with n > 4 vanish in four dimensions. Therefore, the massive gravity 
theory can be written as 

£ = ^ V=9 (fl - . (68) 

where U is given by ( 66 ). 

The theory dehned by ( 68 ) has Minkowski spacetime as a trivial solution, hence 
one can rewrite the metric and the scalars as deviations from hat space, namely 

g^,u = ViMu + 0“ = - 7r“, (69) 

where x°‘ are the usual cartesian coordinates spanning rjay- 

We focus on a convenient limit of Lagrangian ( 68 ) which captures most of the 
dynamics of the helicity -0 mode, but keeps the linear behaviour of the helicity -2 (tensor) 
mode [126]. The limit, called the decoupling limit, is dehned as 

m —)■ 0 , Mpi — )■ cx) , A 3 = m^Mpi = hxed. (70) 

In order to obtain canonically normalized kinetic terms for the helicity 2 and helicity 1 
modes, together with the relevant couplings for the helicity 0 modes, when this limit is 
taken one needs to canonically normalise the helds in the following way 

h^iv —)■ mMpi , n —)■ m^Mpin, 


(71) 
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where we have split the Stiickelberg helds into a scalar component tt and a 
divergenceless vector A^, namely 


+ A,). 


(72) 


The resulting Lagrangian in the decoupling limit is [130] 

£ = - ^ V + A"" (V<L> + x;s + A'W), (73) 

where we have ignored the vector. This decoupling limit theory is a special case of 
the effective theory derived from the Horndeski theory Eq. (52). Thus in this limit the 
theory shares the same properties as the Horndeski theory. 

Cosmology in this massive gravity theory is rather peculiar. If we assume that the 
hducial metric obeys the FRW symmetry only the open FRW solution is allowed 
[132, 133]. We can allow an inhomogeneous hducial metric by choosing the Stiickelberg 
helds of the form [134, 135, 136, 137, 138] 

<P' = 9{t,r) — - (74) 

r 

Then there exists a class of solution characterised by 


a(t, r) = r)r, c, = 1 + 6a. + 12a. ± + 3a3 + 9«1 - 12a, ^ 


3(1 T 303 T 404) 


The Einstein equations are given by 


Gf, = = 


1 T 3q; 3 i 2 q;5 
3(1 + 3^3 ±05)2 


m , 


(76) 


where 


= 1 + 3^3 + 9al - 12a4 • (77) 

The graviton mass plays the role of an ehective cosmological constant. This model 
could realise the idea of technical naturalness. Since the massless limit has more 
symmetry (i.e. GR), the small mass is natural in the sense that the quantum corrections 
only renormalise the mass as —)■ [139]. Unfortunately, however, 

perturbations around these solutions are found to be unstable [140, 141, 142, 143]. 

A natural extension of massive gravity is realised by promoting the hducial metric, 
SO, to a dynamical metric described by another Einstein-Hilbert action. This is known 
as bigravity. The action for the massive gravity model (68) can be extended to the 
bigravity model without introducing a ghost [144]. In this case, it is possible to hud FRW 
solutions for both of the metrics [145, 146, 147]. Unfortunately, again these solutions 
are found to be unstable under linear perturbations [148, 149, 150, 151, 152, 153, 154] 
unless we tune the Planck scale for the second metric to decouple the massive graviton 
from matter [155]. It is also possible to modify the coupling to matter by introducing 
an ehective metric that is made of two metrics. Generally this brings back the BD 
ghost [156]. There is a coupling which removes the BD ghost in the decoupling limit 
[157, 158]. However, linear perturbations around the cosmological solutions are again 
unstable [159, 160, 161, 162]. These massive/bi-gravity theories can be formulated 



Cosmological Tests of Modified Gravity 


23 


elegantly using the tetrad formalism [163]. It has been suggested that there is a matter 
coupling that does not introduce the BD ghost in the full theory although it is not easy 
to hnd a corresponding metric formulation [164] (see also [165, 166]). See Ref. [31] for 
other extensions of the model and a more complete list of references. It remains to be 
seen whether there exist stable solutions where the acceleration is driven by the mass 
and are also distinguishable from ACDM observationally. 

2. 7. Outlook 

Although there has been signihcant progress in developing modihed gravity models and 
many interesting new ideas have been explored to tackle the cosmological constant and 
the late time acceleration problem, we still do not have a consistent theory that is a 
true alternative to ACDM. Our quest for modihed gravity models as an alternative to 
ACDM still continues. 

However, many different models that have been reviewed so far have common 
properties. There are usually three regimes of gravity: 

• There is a length scale above which gravity is modihed. This can be simply 
determined by the cosmological constant, or the Compton wavelength of a massive 
graviton, or the cross-over scale in a braneworld model. 

• Even below the modihcation scale, gravity is still modihed due to the extra scalar 
degree of freedom. In this regime, gravity is described by scalar tensor gravity with 
a 0(1) Brans-Dicke parameter. 

• On small scales there appears a scale below which GR is restored due to the non¬ 
linear interaction of the scalar held. This is often called the screening mechanism. 
This is essential to ensure that the theory passes the stringent Solar System 
constraints. 

This general picture is useful to develop a strategy to construct cosmological 
and astrophysical tests of gravity even though the examples that we have seen suher 
from various theoretical problems [167]. One of the major developments of modihed 
gravity is the screening mechanism, which enables us to modify gravity signihcantly 
on cosmological scales yet to satisfy the Solar System constraints. We will review the 
screening mechanisms in the next section. 

3. Screening mechanism 

3.1. Why do we need a sereening mechanism? 

As we learnt from several examples of modihed gravity models, one of the difficulties in 
explaining the accelerated expansion of the Universe by modifying gravity is that GR 
is tested to high accuracy in the Solar System. In the simplest example of Brans-Dicke 
gravity, once we impose the Solar System constraints, there is essentially no room to 
modify gravity on cosmological scales. One possibility to avoid this problem is to break 
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the equivalence principle. The stringent constraints in the Solar System are obtained 
using objects made of baryons. If the additional degree of freedom only couples to dark 
matter not to baryons, we can evade these stringent constraints while it is possible to 
modify gravity signihcantly on cosmological scales. This is known as interacting dark 
energy models in the Einstein frame where the scalar held is coupled only to dark matter 
| 168 |. 

If we keep the equivalence principle, we need to have a mechanism to suppress the 
modihcation of gravity on small scales. We loosely call this mechanism the screening 
mechanism. The screening mechanism is realised by the fact that the additional degree of 
freedom, which is often represented by a scalar held, obeys a non-linear equation driven 
by the density. The density varies over many orders of magnitude in our Universe. The 
critical density of the universe is given by pcrit = 10“^® g cm“^. The typical density 
inside galaxies is pgai = 10“^^ g cm“^. The density in the Sun is given by psun = 10 g 
cm“^. Thus if we expand the density around the cosmological background, the density 
contrast in the environments where we perform conventional tests of gravity is much 
larger than one. The screening mechanism utilises the non-linearity of the scalar held 
driven by the non-linear density contrast to change the behaviour of the scalar held 
from cosmology to the Solar System. 

3.2. Screening mechanisms 

In this section we follow Ref. [11] and classify the screening mechanisms. A general 
Lagrangian for a scalar held can be written schematically as 

C = dfi, - V{fi) + /3(0)t;, (78) 

where represents derivative self-interactions of the scalar held, V (0) is a potential, 
/3{4>) is a coupling function and Tfi is the trace of the matter energy-momentum tensor. 
In order to avoid the Ostrogradsky ghost associated with higher time derivatives, 
contains up to the second derivative of the held. In the presence of non-relativistic 
matter Tfi = —p, the scalar held’s dynamics depends on the local density of the system. 
Let’s consider the background held 0 which depends on the local density. Around this 
background, the dynamics of huctuations is determined by three parameters: the mass 
m(0), the coupling P{4>) and the kinetic function Screening can be realised 

mainly in three diherent ways utilising these three parameters: 

• Large mass 

If the mass of huctuations w?{(j)) is large in dense environments, the scalar held 
does not propagate above the Compton wavelength and the additional 

force mediated by the scalar held is suppressed. On the other hand, in low density 
environments such as cosmological background, the mass can be light and the scalar 
held mediates the hfth force modifying gravity signihcantly. This idea is realised in 
the chameleon type screening mechanism [44, 45]. 
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• Small coupling 

If the coupling to matter /3{4>) is small in the region of high density, the strength of 
the hfth force generated by the scalar held is weak and modihcations of gravity is 
suppressed. On the other hand, in low density environments, the hfth force strength 
can be of the same order of gravity. This idea is realised in the dilaton [169] and 
symmetron mechanism [170]. 

• Large kinetic term 

If we make the kinetic function large in dense environments, the coupling 

to matter is ehectively suppressed. There are two possibilities to make the kinetic 
term large. One is to assume that the hrst derivative of the held becomes large. 
This idea is realised in the k-mouhage type mechanism [171]. On the other hand, 
in the Vainshtein mechanism [129] , the second derivative of the held becomes large 
in the region of high density. 

In the next section, we look at examples of models that accommodate these screening 
mechanisms. 


3.3. Examples 


3.3.1. Chamleon/Symmetron/Dilaton mechanism This class of model is represented 
by an action (23) in the Einstein frame 





1 

IGttG 







(79) 


The matter helds couple to a metric Due to this coupling, a test particle feels 

the hfth force Vln74((;/)) generated by the scalar held. In these models, the dynamics of 
the scalar held is determined by the local density dependent ehective potential 


Veff = V{fi) - [A{4>) - l]T/. (80) 

The dynamics of the scalar held is characterised by the mass of the scalar held aronnd 
the minimum of the potential (j) = and the conpling fnnction 


m 


= VM) /3 = Mpi 


d\n.A 


d(j) 


(81) 


Below we show typical choices of the potential and the coupling function to realise 
the screening mechanisms: 

71(0) = 1 + V(0) = —— chamaleon, (82) 

Mpi (p^ 

M/) = ^ dilaton, (83) 

M/) = 1 + ^ symmetron, (84) 

where the mass scale M is a parameter of these models. In the chameleon mechanism, 
the mass becomes large in the region of high densities while in the dilaton and 
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symmetron models, the coupling fi becomes small in high dense environments. Note 
that the second derivative of the effective potential does not depend on the enegy density 
explicitly in chameleon models as A{(j)) is lienar in fi. The density dependence of the 
mass comes from the fact that the minimum of the potential 0 depends on the density. 

In all these models, if we consider an object with the gravitational potential 
the screening of the hfth force happens when the gravitational potential |Tn| exceeds 
some critical valne iTjvl > tjjc- We will show this in the case of the chameleon mechanism 
below. 

3.3.2. K-moflouge/Vainshtein mechanism The other class of the model utilises the 
non-linear derivative interactions to screen the hfth force. Typical choices of the 
fnnctions in the Horndeski action (46) 

^ = j d^x^ ^/? + ir(0,X)-G3(0,W)n0 +SU9,u). (85) 

to realise the screening mechanisms are given below: 

X) = X + , k-mouhage, (86) 

K{(j),X) = X, G'3(0,X) = ^X Vainshtein. (87) 

The screening condition for a spherically symmetric object with M and radius R 
is determined by the solntion for the scalar held fi. Again let us consider an object 
with the gravitational potential |4/Ar| = GM/R. In the k-mouhage mechanism, the 
screening of the hfth force happens when the hrst derivative of the potential, i.e. the 
gravitational acceleration exceeds some critical value |(9TAr| > Ac. On the other hand, 
in the Vainshtein mechanism the screening operates when the second derivative of the 
potential, i.e. the spatial curvatnre, exceeds some critical value |i9^TAr| > Iff. We 
will show this explicitly for the Vainshtein mechanism below by solving the scalar held 
equation. 


3 . 4 . Chameleon screening 

Using a simple example, we show how the chameleon mechanism snppresses the scalar 
held [45]. Let’s consider a scalar held with mass m that is coupled to matter. We assume 
that the function A(0) can be approximated as ^(0) = 1 -|- fi(j)/M-p\. In the case of BD 
gravity 


2(3 + 2n;BD)’ ^ ^ 

which gives = 1/6 for f{R) gravity. In the Einstein frame, the geodesic eqnation of a 
test particle is modihed due to the coupling between the scalar held and matter. From 
the geodesic eqnation, the acceleration of the particle is obtained as 

a = -VT - = -VT - /3V0, (89) 

d<p 


-V0 = -VT - /3V0, 
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and the equation of motion for the scalar field is given by 

= SirCfip, 


(90) 


where we assume does not depend on time. A spherically symmetric solution outside 
the matter source is given by 


0(r) = 


2fiGM 


(91) 


where M is the mass of the object. For scales smaller than the Compton wavelength, 
we can ignore the Yukawa suppression in Eq. (91). The acceleration of the particle can 
be obtained by substituting Eq. (91) into Eq. (89) as a = —Ge^M/r"^ where 

'A + 2ujbd\ 

.3 + 2ojbd ' 


Gen = G 1 + 2/?q = G 


(92) 


Thus the test particle feels the modified Newton constant as in BD gravity without 
potential (see Eqs. (17) and (18)). On the other hand, due to the Yukawa suppression 
at length scales larger than the Compton wavelength m~^, the scalar field does not 
propagate and the test particle feels normal gravity. The idea of chameleon screening is 
to make the mass dependent on the density so that the scalar field does not propagate 
in the dense environment. 

To realise this idea, we need to introduce a potential to the scalar field. The 
equation of motion is generalised to 


VV = V;'|,(»J) = + 8!rG/3p. 


(93) 


By choosing the potential Y(0) appropriately, we can realise the situation where the 
mass of the scalar field is large for larger densities. 


log(l) log(l) 



Figure 3. The effective potential —VeS inside and outside the matter source. We 
can understand the spherically symmetric solution as a dynamically time dependent 
field governed by — Kff- The red curve is the bare potential —V and the blue curve 
represents the density dependent contribution. The back curve is the total effective 
potential —Ves. 

Let’s consider a spherically symmetric matter source with a constant density within 
the radius r = R. Inside the source, the scalar field is trapped at the minimum of the 
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effective potential = 0, 0 = fis- Outside the source, we assume that we have the 
background density poo ^ Ps and the scalar held is again trapped by the minimum 
0 = 0OO with mass moo- The dynamics of the scalar held is governed by the equation 

0(r)" + -0'(r) = O'(0) + 8nG/3p. (94) 

r 

By identifying r as time t, this is equivalent to the time evolution of the scalar held with 
a potential —Ves (Fig. 3). The scalar held starts from the “maximum” of the potential 
(j) = fis and then rolls to another “maximum” 0 = fioo- The solution for the scalar held 
outside the matter source is obtained as 0(r) = — (C'/r)e“™'°°” + fioo- From the radius 
r = R until the scalar held settles down to 0 = 0^ at radius r = Vscr, the scalar held 
dynamics is driven by the change of the density so we can ignore the contribution from 
the potential. Then the equation of motion for the scalar in this regime is given by 

0(r)" + -0'(r) = SttGPp, (95) 

r 

which admits the solution 0 = {apc/8Mpi)r^/2 + A/r + B. The integration constants 
A, B and G, and the radius rs^r are determined by matching the held and its hrst 
derivative at r = R and r = Tscv We obtain the solution outside the source as 


0(r) 




(96) 


where 

^R _ 0OO 0s dd Tscr 

^ 60Mpi|T^| R ’ ^ ^ 

and is the gravitational potential of the object |Tjv| = GM/R. If 3AR/R > 1 this 
solution is not valid and the solution becomes Eq. (96) with 3AR/R = 1. The condition 
AR/R 1 is called the thin shell condition. If this condition is satished only the mass 
within the thin-shell of size AR contributes to the force because in the interior of the 


source the scalar held is heavy and the scalar force is Yukawa suppressed. Note that we 
naively say that the scalar held is screened in dense environments. This is in fact not a 
very accurate statement. What matters is the gravitational potential of the object not 
the density. See Fig. 4 for a schematic picture of the thin shell solution. 

Although we derived the solutions in the Einstein frame where the dynamics of 
the scalar held is easier to understand, it is straightforward to translate the solutions 
into those found in the original (Jordan) frame where the scalar held is not coupled to 
matter directly while it couples to gravity non-minimally (see Eq. (12)). The solutions 
for the metric and scalar held perturbation are given by 

= dvrGpp, T = 77 “^$, (98) 


p = l + 2fil 


1 2 _ o2^^di 

1 , oo2 ’ PeS P TO ’ 


when the thin shell condition is satished. The scalar held in the Jordan frame 0 is 


related to 0 as 0 = —200/Mpi for 0 1. In particular, for f{R) gravity we have, 

i/j = fR = — v^2/3(0/Mpi) (see Eq. (26)) and the thin shell condition can be written as 

I/tjoo —/i?s| < 2 |TAr|. (100) 
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Figure 4. A schematic picture of the thin shell profile of the scalar field. 


Let’s consider the Solar System constraint [46]. In the Sun, the density is p = 10 
g cm“^ while the density of the Milky way galaxy is typically g cm“^. Thus the 

scalar held in the Solar System is well suppressed compared with the one in the Milky 
way galaxy fisoiar ^ figni- If the thin shell condition is satished, is small so we can 
approximate p as p — 1 = The constraints |p — 1| = (2.1 ± 2.3) x 10“^ gives a 

constraint on the galaxy held 


^ < 10 -“, ( 101 ) 

iWpi 

where we used the fact that the Milky Way’s potential is approximately given by 
HIAlgal ~ 10“®. This is a model independent constraint but in order to translate this 
constraint to the mode parameters, we need to interpolate this constraint to the scalar 
held in the cosmological background ficos- For this, we need to specify the potential. We 
consider a potential of the form 


V = A- 


M, 


PI 


1+n 


( 102 ) 


as in the f{R) gravity given by Eq. (24) (see Eq. (25)). We assume that the scalar held 
is at the minimum of the potential 

-(l+n) 

(103) 


M, 


PI 


(n + l)/3p Y 
nM^ ) 


Given that the ratio of the cosmological background and the Milky way galaxy is 
Pgal/ Pcos 10 ®, the ratio of the scalar held is given by 

-(l+n) 

= (104) 


pc 


Pgal 


Using this relation between 0gai and ficosi we obtain the constraint on the cosmological 
held (pcos 

To derive the constraint, we assumed that the galaxy held is at the minimum of the 
potential. In fact, if the Milky way is not screened the dynamics of stars are disrupted. 
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Imposing the thin shell condition for the Milky Way galaxy, we obtain the constraint 
fificos!dM-p\ < 10“® where we again have used "^Ngai ~ 10“® and 0gai ficos- This 
constraint is much more stringent than Eq. (101) for n > 0. Note that to derive the 
constraint that the Milky Way galaxy is screened, we assumed that the Milky way is an 
isolated system in the cosmological background. There is a possibility that the Milky 
Way galaxy is screened by the Local Group with a potential \'^n\ ~ 10“^. In this case 
the constraint is relaxed to be fifioo/Mpi ~ 10“^. 

For the potential given by (102), laboratory tests give much weaker constraints 
compared with the cosmological constraints [13]. See [13, 11] for laboratory tests of the 
chameleon models. 

The constraint that the Milky Way with the potential [TtvI = 10“® is screened 
has important implications. Firstly, this unfortunately excludes the possibility of self¬ 
acceleration [172]. The change of 0 between the cosmological density and the galactic 
density Pgai = lO^pcos is suppressed by the thin shell condition /3(0cos —0gai)/Tfpi < 10“®. 
If we translate this to the cosmological evolution of the scalar held, this implies that 
the scalar held barely moves between 0 < < 1 where is the redshift (see [173] for 

caveats). Thus the scalar held cannot be responsible for the acceleration. This also 
excludes the possibility to modify gravity on cosmological scales. The mass of the scalar 
held around the minimum in the cosmological background is given by [172, 174] 


Pc 


Mpi 


Ml 


PI 


< 10-®M 


(105) 


Thus there is no modihcation of gravity on scales larger than 1 h ^Mpc. Another 


IT 


N\ 


< 


implication is that GR is modihed only with objects with a shallower potential 
10“®. We will discuss how we test models with the Ghameleon mechanism using these 
objects later. 

Finally we discuss briehy about the quantum corrections. The one-loop Goleman- 
Weinberg correction to the potential grows as Imposing the condition that the 

quantum corrections are under control at laboratory density piab ~ 10 g cm“^, the 
constraint on the mass of the chameleon is given by [175] 


< 


Mpi J 


1/6 


= 7.3 X 10-3 


fipiab 


1/3 


eV. 


(106) 


10 g cm“3^ 

This upper bound is already in tension with the laboratory constraints on the hfth force. 
This also means that the cut-off scale of this theory is rather low Apv ~ lO-^eV ~ ■ 

The coupling to matter also generates corrections to the potential but this contribution 
is suppressed by the cut-off scale of the theory Apv as Am/m = so this is 

totally negligible. 


3.5. Vainshtein mechanism 

We consider the simplest cubic galileon model to understand how the Vainshtein 
mechanism operates. We use the equation of motion for the scalar held in the DGP 
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model, Eq. (30): 

3/?VV + ’■d(VV)^ - (ViV,i#.)(V‘V^i^)l = SnGSp. 


(107) 


For the spherically symmetric case, it is possible to integrate the equation once. Outside 
the matter source, the solution obtained is [52] 


d(j) 

dr 



rv = 


Sr^rg 

9/32 


1/3 


(108) 


where = 2GM. Below the Vainshtein radius ry, the non-linear interaction suppresses 
the hfth force and we recover GR solutions with small corrections [51, 50]. The solutions 
for the metric perturbations can be obtained by subsituting the scalar held solution (108) 
in the limit r <C ry into Eqs. (14) and (16) as 


$ 


GM ^ 
r 






(109) 


where the upper (lower) sign corresponds to a positive (negative) /3, i.e the normal 
(self-accelerating) branch. On the other hand, outside the Vainshtein radius, the linear 
solution is recovered [51, 50, 176] 


$ = 






( 110 ) 


See Fig. 5 for a schematic picture of the Vainshtein solution. Note that the condition to 



Figure 5. A schematic picture of the Vainshtein solution (107). The Vainshtein radius 
ry at which the non-linear interaction of the scalar field becomes important is much 
larger than the gravitational radius rg at which the spin-2 graviton becomes non-linear. 
The cut-off scale of the theory Ag ^ is typically very large. 

realise the Vainshtein mechanism r ry can be written as > r“^, conhrming 

that the Vainstein mechanism operates when the spatial curvature exceeds a critical 
value. 

The Vainshtein radius ry is much larger than the gravitational length for Vc Ho. 
For the Sun, r^ = 2.95km while the Vainshtein radius is ry = 130 pc. Thus inside the 
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Solar System, the solutions for the metric are well approximated by Eqs. (109). The 
perihelion precession per orbit Afi is calculated as [51] 

1/2 

( 111 ) 


A , STTro Stt 

A0 = 27r +-^ ± — 

r 2 


2r‘^rg 


where the second term is the usual Einstein precession and the last term is due to the 
hfth force. The angle of perihelion advance due to the hfth force during one orbital 
period is then given by 

1 

arcseccond per century. (112) 


Ad) 3 
— = — = 5 X 10 
P SVr 


-4 


5Gpc 


The interesting feature of this prediction is that the precession rate is independent 
of mass so it is universal. Combining the observed precession of planets in the Solar 
System, the constraint was obtained as Afi/P < 0.02 arcsecond per century, which 
gives the constraint Vc > 130 Mpc [177]. Another constraint comes from Lunar Laser 
Ranging experiments. These experiments put constraints on the correction to the 
Newton potential < 2.4 x 10“^^. Using Eq. (109), Vg = 0.886 cm for the Earth 

and r = 3.84 x 10^° cm for the Earth-Moon distance we obtain Vc > 162 Mpc [178]. 
This constraint will be improved by a factor of 10 by the Apache Point Observatory 
Lunar Laser-ranging Operation (APOLLO). See Ref. [179] for a review. These results 
are obtained by treating the planets or the Moon as a test body. This assumption is 
not necessarily valid. For example for the Earth-Moon system, it was shown that there 
is a correction to the universal precession rate (112) due to the non-superimposability 
of the held that depends on the mass ratio of the two bodies, (although the effect is 
small (4%) for the Earth-Moon system) [180]. See Ref. [181] for the constraints from 
laboratory tests. 

The most general action describing the Vainshtein mechanism can be obtained from 
the Horndeski action as we described in section 2.5. From the effective action (52), it 
is easy to hnd static spherically symmetric solutions [90, 89]. We consider the following 
conhgnration 


flu = -2$, hij = 0 = 00 + 7r(r). (113) 

The held eqnations yield 

P{x, A) = fA(r) + (I + + 6 q;,^ - 3fiA(r)'jx^ 

+ (^z/ + 2a2+ 40e)x^-30V = 0, (114) 

y(r) = fix^ + A(r), (115) 


where we dehne 

x{r) 

A{r) 


IP , . Mpi$' 

Af7’ “ 7f7 

1 M{<r) 

MpiAg 


Mpi T' 
Ag r ’ 


(116) 


(117) 
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The function M(< r) represents the mass of a spherically symmetric, pressure-less 
matter source, up to a radius r. Outside the surface of the matter source, A{r) can 
be written as 



/ M \ 1/3 
V SttMpi Ao ) ’ 


( 118 ) 


where ry is the Vainshtein radius depending on the total mass M of the source. 

Stability of these solutions can be studied by considering small perturbations around 
the background tt = 7ro(r) -|- ip{t, r, O), <h = <h(r) -|- r, O), T = 'h(r) -|- (5T(f, r, O) 
where O represents angular coordinates. By expanding the effective action up to the 
second order in these small perturbations we hnd [90] 


where 



d^x 




Kr{r){drip) 


KQ{r){dnipf , 


(119) 


Kfir) = 2dxP{x,r)\x=xo, 


( 120 ) 


Kfir) = 


3^2 + 6Q;,^)a; + 12 q;A 

-|-18(z/ -|- 2 q;^ -|- 4:^)x‘^ + 12(10q;/3 + w)x^ + 


Koir) = 


r^|?7 -I- 6.^^ -I- 4(yn -|- 6a^)x 


2 r dr 

+ 2a^ + 4^)x‘^ — 


( 121 ) 

( 122 ) 


In addition to the quadratic action for scalar perturbations, there is also a coupling 
between metric perturbations and the scalar perturbations. The stability of fluctuations 
depend strongly on the presence of two particnlar couplings in the Lagrangian (52), one 
given by and the other by , which have a and fi as coefficients, 

respectively. We will consider their conseqnences in what follows, by analysing each 
relevant case separately [90]. 


• a = (3 = t). 

In this case the action (52) reduces to the galileon action introdnced by Ref. [54]. 
Interestingly, it was fonnd that the stability conditions force the the speed of sonnd 
for fluctnations propagating radially to be greater than one. On the other hand, 
if the fourth order galileon term is included, then Kq Kt,Kr and the 
propagation of angnlar fluctuations is extremely sublnminal, invalidating a quasi¬ 
static approximation. 

• CK 7 ^ 0 and (3=0. 

A new feature of this class of models is the disformal coupling between matter 
and the scalar, namely d^ndiyiiT^’^. As pointed ont in Ref. [182], this conpling has 
deep implications for the stability. Inside a matter sonrce, the coupling introdnces 
a kinetic term proportional to ap{dt(pfi. If a < 0, the fluctnations behave as 
ghosts, forcing us to choose a > 0. In the case of massive gravity, there is only one 
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dimensionless parameter, a = 1 + 80:3 (in the theory where (3 = + 4 q ;4 = 0 ), 

leading to the so-called restricted galileon [182], In this case, there is no Vainshtein 
solution that can be connected to the asymptotically flat solution. Around the 
non-flat asymptotic solution, the radial propagation is subluminal, even though the 
extreme subluminality of the angular propagation persists [182], 


• a 0, /3 ^ 0 

In this case, the coupling between the metric and the scalar gives a different 

picture from the previous cases. Inside the Vainshtein radius the solution that 
reduces to GR near a matter source only exists if > 0, and the solution is given 
by [183, 184] 


Xo = ± 



(123) 


Since x is constant and A{r) 3> 1 inside the Vainshtein radius, all functions Ki 
in (120) are dominated by the ?/ ~ A(r) 3> 1 contribution in the limit r ry, 
as shown from its definition (116). In the limit A{r) 3> 1, we can also ignore 
the coupling between the scalar and metric perturbations. By taking into account 
only the contributions depending on A = {ry/rfi outside the surface of the source 
r > Tg, we And in this limit 


Kt = 0, Kr = -l2/3A{r)xo, Kn = 6l3A{r)xo. (124) 


Thus inside the Vainshtein radius r ry, but outside the source surface, the speed 
of the fluctuations are always superluminal, for both radial and angular directions. 
Moreover, given the fact that Kq and have opposite sign, all solutions are 
unstable. 


These results indicate that it is not always possible to have a successful Vainshtein 
mechanism. For example, there is no stable Vainshtein solution that connects to the 
asymptotically flat spacetime in massive gravity theory discussed in section 2.6. See 
Refs [185, 186, 187, 187, 188, 189, 190] for the studies of the Vainshtein mechanism in 
the Horndeski, beyond Horndeski and bigraviy theories. 

Away from the decoupling limit, numerical methods are required to examine the 
solutions. Refs. [191, 192] found solutions featuring the Vainshtein mechanism in the 
complete theory in the framework of a non-linear extension of Fierz-Pauli massive 
gravity. See [193] and references therein. There have been several studies of the 
Vainshtein mechanism away from the spherically symmetric case (see for example 
[194, 195]) and weak gravity [196, 197]. 

There are also various theoretical problems associated with the non-linear 
interaction terms. The first problem is the superluminality. As we discussed earlier, 
if we consider perturbations around spherically symmetric solutions, the radial speed 
of propagation is larger than the speed of light in galileon models. This is a generic 
feature of galileons as the derivative interactions change the structure of light-cone for 
perturbations. The consequence of this superluminality is still debated [198, 199, 200] 
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(see [31] for a summary of discussions). Another problem is the strong coupling problem. 
The quantum corrections become important at energy scales larger than A 3 . In the 
DGP case, A 3 ^ = (Mpir- 2 )-V 3 _ 1000km to satisfy the Solar System constraint. This 
strong coupling scale depends on the background as seen from Eq. (119). Around the 
Vainshtein solution, the strong coupling scale is renormalised [72] 


Af =(yl 


ry \ 3/4 




Aq 


(125) 


which reduces the strong coupling scale to 1cm around the Earth. This low cut-off 
scale is a generic problem in theories with galileon interactions such as massive gravity 
[201]. Although quantum corrections are not under control above A 3 , the galileon terms 
themselves do not get renormalised upon loop corrections so that their classical values 
can be trusted quantum mechanically [ 202 ]. Galileons also enjoy self-duality, i.e. various 
galileon terms are remapped to other galileon terms by a held dependent coordinate 
transformation and a held redehnition [203]. It is even possible to hnd galileon terms 
that can be mapped to a free theory. Galileon duality was used to argue that the 
strong coupling and the superluminality problem are not fundamental [204] while it was 
also argued that the Ultra Violet (UV) completion inhuences the low energy theory 
signihcantly [205]. Glearly more work seems to be warranted to understand the UV 
properties of galileon models. 


4. Cosmological tests of gravity 
4 -1. Consistency test 

Armed with the theoretical knowledge of modihed gravity models, we now discuss how 
we test gravity on cosmological scales. The hrst approach is to test specihc modihed 
gravity models using observations. The models discussed in section 2 have been tested 
intensively by observations (see [43] for f{R) gravity, [30] for the DGP models and 
[206, 207, 208] for covariant galileons). We will not attempt to review observational 
tests of individual models. 

At the background level, the Friedman equation in modihed gravity models can 
always be recast into the Friedman equation in GR with dark energy 

= + (126) 

By tuning the equation of state for dark energy w^e = Pde/ Pde h is always possible to 
mimic the background expansion of the Universe in modihed gravity models by dark 
energy. The equation of state is commonly parametrised as [209] 

WAe{a) =Wo + Wa{l-a). (127) 

This degeneracy is broken if we include observables determined by the structure 
formation assuming that dark energy is smooth. The most commonly used 
parametrisation of the structure growth is to parametrise the growth rate, the 
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logarithmic derivative of the dark matter density perturbation with respect to the scale 
factor, [ 210 ] 


/ = 


d log 6 „ 




(128) 


dloga 

where 6 m is the density fluctuation of dark matter. In ACDM, 7 is approximately given 
by 7 ~ 0.55. In dark energy models 7 is completely determined by the equation of 
state Wde- However, in modihed gravity models, 7 is not related to the equation of 
state. Thus if our universe were described by a modihed gravity model but we tried 
to ht various observations including those measuring the background expansion of the 
Universe and the growth of the structure, we would hnd inconsistencies [211, 210]. Fig. 6 
demonstrates this fact [ 212 ]. 




Figure 6. Left: Forecasted 1 cr constraints on the dark energy equation of state (127) 
for the Dark Energy Survey [4]. From the largest to the smallest ellipse, the probes 
considered here are Baryon Acoustic Oscillations (black), Supernovae (green), cluster 
counts (magenta) and weak leasing (blue). The true model is assumed to be ACDM. 
Right: Same as the left hgure but the true model is assumed to be a modified gravity 
model with 7 = 0.68. Due to the incorrect attempt to ht a GR + dark energy model 
to the data, the centres of the weak leasing and cluster counts ellipses have moved. 
From [212]. 


We can formulate the consistency test of ACDM formally using linear cosmological 

perturbation theory [213, 214]. In ACDM the basic equations that govern the expansion 

of the universe and the growth of structure are given by 

SttG , . 

H = pt, (129) 


—4/ = AtiG pt6t 1 (130) 

a"* 

where pt is the total energy density and 6 t is its fluctuation in the comoving gauge. Here 
we are working in the Fourier space for perturbations. These equations are supplemented 
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by the conservation of energy momentum tensor. Let’s assume that the dominant 
component of the energy density in the late-time universe is dark matter, px = Pm, 
dr = dm- The continuity and Euler equations give 

L + -Om - 3$ = 0, (131) 

a 

u2 

dm + H9m -T = 0, (132) 

a 

where 9m is the velocity divergence 9m = '^ ■ ^m- Under horizon scales, we can neglect 
<h in the continuity equation and we hnd a second order equation describing the growth 
of dm 

dm + 2Hdm = -’^'^. (133) 

This equation describes how the density perturbation grows due to the effect of the 
Newtonian potential T. Thus the modihcation of gravity predicts a different growth 
rate. 


There are a number of cosmological observations that measure different geometrical 
quantities (see [11, 15, 14] for reviews). The Cosmic Microwave Background (CMB), 
Supernovae (SNe) and the Baryon acoustic oscillations (BAO) measure the distances 
in the background universe, hence the expansion history of the Universe H{z). Weak 
gravitational Tensing (WL) measures the distortion of galaxies due to the deflection of 
light, which is determined by the lensing potential $ -|- T. The time variation of the 
lensing potential changes the temperature of the CMB photons through the Integrated 
Sach-Wolfe (ISW) effects. Peculiar velocities of galaxies follow the velocity of dark 
matter 6 m, which changes the distribution of galaxies in the redshift space. This Redshift 
Space Distortion (RSD) can be used to measure the velocity of dark matter. The galaxy 
distribution is a biased tracer of the underlying dark matter distribution. Although it 
is very difficult, in principle, we can reconstruct the dark matter density perturbation 
dm from the distribution of galaxies. 

By combining the Friedman equation and the Poission equation, we obtain the 
following equation [214, 215] 

2k^ ($ + T) - T 


a{k, t) = 




dr, 


= 1 . 


(134) 


This is the identity in ACDM. Let’s examine each term in this expression. First {k/a) is 
the physical wave number of the perturbations that we are interested in. The background 
expansion H can be measured from CMB, SNe and BAO. We split <h into T -|- <l> and 
— T. 4/ -|- <h is the lensing potential which can be measured from WL and ISW. On the 
other hand, the Newton potential 4/ governs the dynamics of galaxies. The peculiar 
velocities of the galaxies are determined by the Newtonian potential /c^4/ = d{a9m)/dt 
(Eq. (132)) thus using RSD, we can reconstruct T. Finally, from the galaxy distribution 
dg we can in principle reconstruct dm if we know bias b = dg/dm- This means that we 
have just enough observables to check the identity (134). This a very strong consistency 
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condition of ACDM: a{k, t) is unity at any time and space in ACDM. Any modifications 
to ACDM lead to deviations of a{k,t) from unity. 

The problem with this consistency relation is that it requires the measurement of 
the dark matter density perturbation 6m- This is difficult due to the uncertainties in 
galaxy bias. A less ambitious test is to combine weak lensing and peculiar velocity 
measurements only. The Eg parameter is introduced as follows [216] 

V2($ + T) 


Eg = 


-3H^a-^e„ 


(135) 


In ACDM, Eg = klm/f- This estimator can be constructed directly from observations 
for example using the lensing-galaxy cross power spectrum and the velocity-galaxy power 
spectrum [217]. 


4-2. Parametrisations 

The degeneracy between dark energy models and modihed gravity models exist in 
principle for perturbations. One can find GR models with dark energy that have 
anisotropic stress and variable sound speed, which can in principle mimic modihed 
gravity models [218]. Mathematically, this is always the case as we can dehne an effective 
energy momentum tensor to absorb any effects of modihcation of gravity 

Ggu = 87rG(T^i, -|- Eg^i^). (136) 

At the background it is enough to specify the equation of state we = Pe/pe for the 
huid whose energy-momentum tensor is given by Eg^y. For linear perturbations, we need 
to specify the sound speed c| and the anisotropic stress tie- Once we hnd deviations 
from ACDM, we need theoretical models which predict these unknown quantities to 
distinguish between dark energy models in GR and modihed gravity models. For 
example, modihed gravity models generally predict a large anisotropic stress tie ~ Pe6e, 
which is difficult to hnd in physical huid models. Thus it is important to go beyond the 
consistency test of ACDM and hnd precisely how the deviations from ACDM appear. 

In general, deviations from ACDM can be parametrised by two functions of 
space (wavenumber k) and time (redshift z) that characterise the relation between 
the geometry and matter. For linear perturbations, one can parametrise the relation 
between the Newton potential to the density perturbations and the relation between 
the two metric perturbations [219, 220] 

k‘^^> = AnGp{k^z)ofipm6m-i (137) 

I = >)(*;. z)- (138) 

Another useful parametrisation which can be constructed from n and rj is to parametrise 
the relation between the lensing potential and the density perturbation 

^2(^ -|- $) = 3nG'P{k^z)a?Pm6m- 
where S = /i(l -|- rj)/2. 


( 139 ) 
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These functions are unity fi{k,z) = rj{k,z) = T{k,z) = 1 in ACDM. If we modify 
gravity they deviate from one. For example, in the BD gravity with a mass term, these 
functions are given by 


S(fc, a) = 1, /i(fc, a) = 


2(2 + ubd) + rrraf /k 


(140) 


2>+ 2ujbd +/k"^ ’ 
under horizon scales. Even in dark energy models, if dark energy clusters but if we 
do not know their existence, this looks like a modihcation of gravity for dark matter. 
Assuming that dark energy has no anisotropic stress, we predict 

PdeSde 


ja{k, a) = T{k, a) = 1 + 


Pn 


(141) 


where is the density perturbation of dark energy. In this way, theories beyond the 
standard ACDM model predict distinct paths in the {fi, S) plane [221]. 

The question is how well we can constrain these two functions of time and space 
from observations. There are many attempts to put constraints on /i and S (or rj) 
[222, 223, 224, 225, 226, 227, 228, 229]. Given the accuracy of current observations, we 
are forced to assume particular time and space dependences of these functions. One 
parametrisation is to assume that these functions are scale invariant and their time 
dependence is determined by the density parameter associated with the cosmological 
constant [228] 


fi{a, k) = /ioDL(a), S(a, k) = SoDL(a), (142) 

The Planck 2015 paper [230] gave constraints on /tq and Sq as shown in Fig. 7 by 
combining Planck CMB observations with various external data (BAO, WL, RSD). The 



Figure 7. 68 % and 95% contour plots for the two parameters Eq and /tq obtained 
from Planck CMB measurements combined with WL, BAO and RSD. From [230]. 

ACDM model is in tension with the data at 3a level if the Planck data is combined with 
WL and RSD/BAO. This tension is reduced to 1.7a level if we include the CMB lensing. 
The WL and RSD measurements will be dramatically improved in the next hve years. 
It is very interesting to see whether this tension remains with improved measurements. 
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In order to improve the constraints, it is important to combine RSD and WL. RSD 
is determined by the peculiar velocity of galaxies thus it is sensitive only to fi{k, z). On 
the other hand, the weak lensing is determined by the lensing potential so it is sensitive 
to T,{k,z) and also jaik^z) as the growth of the dark matter density is determined by 
/i(/c, z). Thus combining these two probes it is possible to break the degeneracy. Fig. 8 
demonstrates this point where the time dependence of /i(a) and S(a) is assumed to be 
/i(a) = /rs(l + a^),S(a) = + a^) (they are assumed to be scale invariant) [225]. 

This hgure also demonstrates that the accuracy of the constraint strongly depends on 
the assumed time dependence. 
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Figure 8. Constraints on /r(a) and S(a). In the left (right) panel, the time dependence 
is assumed to be s = 1 (s = 3). The constraints are obtained using WL (CFHTLS) 
and RSD (SDSS DR7 LRG). From [225]. 

If we do not wish to assume any theory, i.e. the functional forms of these functions, 
the best strategy is to simply make bins in terms of redshifts and wave-numbers and treat 
fi and T] as free parameters in each bin. The problem is that errors on these functions are 
highly correlated. Principal Component Analysis (PCA) provides a way to de-correlate 
the errors by creating linear combinations of these parameters [231, 232, 233]. Let’s 
consider m z-bins and n k-bins. There are 2 x m x n parameters for p and 7 . The 
covariant matrix is given by 

Cij = {{p - Pi){p - Pj)), (143) 

where pi is the best £t values. The covariant matrix will be non-diagonal as the 
individual pixels of p and p are highly correlated. We can diagonalise the covariant 
matrix 

C = IT^AIT, Aij = XAj. (144) 

Using the matrix W, we can construct the linear combinations of pp 
“i = - ft)- 


( 145 ) 
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The error on at is given by A* and they are nncorrelated. The parameter pi is formally 
expanded by the eigenvectors Wj{k,z) 

Pi{k,z) = pi + '^ajWj{k,z). (146) 

i 

Fig. 9 shows the nnmber of PCA eigenmodes and their uncertainties for an Euclid- 
like survey, which combines imaging (2D, WL) and spectroscopic (3D, RSD) surveys 
[234], We can clearly see that by adding the information from RSD to WL, the 
constraints on the eigenmodes are significantly improved. 
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Figure 9. The number of PCA eigenmodes and their uncertainties for an Euclid-like 
survey. Here 7 = 77 . The lower panel shows the improvements of the errors when 3D 
(RSD) information is added to 2D (WL) measurements. From [234] (published on 19 
August 2013 ©SISSA Medialab Sri. Reproduced by permission of lOP Publishing. 
All rights reserved). 


4-3. Effective field theory approach 

Although the parameterisations Eq. (138) are the most general form for the linear 
perturbations, there is still a gap between the constraints on p and S and the predictions 
of theoretical models. Firstly, these are the parameterisations of the solutions for metric 
perturbations and the density perturbation. In principle we need to solve a complete 
set of equations for linear cosmological perturbations in a given theory to find solutions 
for p and S. Thus these functions depend not only on the parameters of theories but 
also initial conditions. Secondly, in general, the modified linearised Einstein equations 
do not have the form given by Eq. (138) except for the case in which the quasi-static 
approximation holds. 
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Several approaches have been developed to make the connection to theories more 
transparent. One approach is to parametrise the linearised modihed Einstein equations 
directly [235, 236, 237, 238, 239, 240, 241, 242], Another approach is to build the 
effective action [243, 244, 245, 246, 247, 248, 249, 250, 251] following the effective field 
theory approach developed for inflationary models [252, 253]. Here we follow the latter 
approach. We consider a single scalar field 0(f) on the FRW background. The presence 
of the scalar field breaks the time diffeomorphism. We can use 0(f) to define a preferred 
time slicing 0 =const. We construct the action using the invariant quantities under the 
spatial diffeomorphisms x* —)■ + It is natural to use the intrinsic curvature and 
the extrinsic curvature of the 0 =const. hypersurface. The extrinsic curvature is 
defined as K^i, = aUy where is the unit vector perpendicular to the hypersurface 




d„>l> 


(147) 


and — Uutiii is the projection tensor. In addition, is invariant under spatial 

diffeomorphisms. If we demand that the equations of motion are at most second order, 
the most general action up to quadratic order can be written as [245, 246] 

M2{tY 


S^Jd^ 

mfitf 


,-ffm 


A{t) - c{t)g^^ + 


-{Sg 


00\2 


6K6g^^ - mi{t){6K^ - + 


fhfitfi 


^^^R5g 


00 


(148) 


where /(f), A(f), c(f), M 2 (t), msit), mfit) and fhfit) are free functions of time. From this 
action, it is possible to derive the modified linearised Einstein equation. These equations 
have been implemented to the linear Einstein-Boltzmann code [254, 255, 256]. 

Based on the effective field theory, it was found that there exists a large model 
space that mimics ACDM on all linear quasistatic subhorizon scales as well as in the 
background evolution [257]. The effective field theory approach allows us to go beyond 
the quasistatic approximations and test these models using the measurement of the 
relativistic contributions to galaxy clustering [258]. 

It is clear that this approach has a strong connection to the ADM formulation of 
the Horndeski theory discussed in section 2.5. In fact in the Horndeski theory, there 
is a special relation between and fh^, = fh^. This relation is broken by beyond 
Horndeski theories. In fact this was the way beyond Horndeski theories were discovered 
|82]. 


Under the quasi-static approximations, it is possible to derive the form of /i and g 
in the case of = rfii [259, 260, 261] 

'pi{a)+p 2 {a)k‘^\ ( l+pfia)}C 


g{k,a) = 


1 +P3{a)k‘^ 


p{k, a) = 


pfia) +pYa)k‘^ 


(149) 


where pfia) are free functions of time determined by the six functions in the effective 
action. See Refs. [262, 262] for related approaches. 

Although the theoretical prior such as the requirement that equations of motion are 
at most second order significantly simplified the form of the two functions fi{k, z) and 



Cosmological Tests of Modified Gravity 


43 


r]{k, z), there are still five free functions of time. The PCA analysis has been performed 
for these functions. Due to the degeneracies, it was shown that these functions are not 
constrained individually even by a future weak leasing survey. On the other hand, about 
10 eigenmodes of the combination of all pa will be measured with uncertainty smaller 
than 1% of the prior by the Large Synoptic Survey Telescope (LSST) project [263]. One 
approach is to reconstruct p,{k,a) and 'y{k,a) using eigenmodes of pfia) following the 
approach to reconstruct the equation of state w{a) using the smoothness prior on pfia) 
functions [264, 265]. 


4-4- Quasi non-linear scales 


On linear scales, it is possible to develop model independent approaches as we described 
in 4.3. However, once the non-linearity becomes important, it becomes difficult to 
develop model independent parameterisations of the non-linear eqnations. On qnasi 
non-linear scales, it is still possible to develop a general framework based on pertnrbation 
theory [266]. In models with a screening mechanism, the scalar field eqnation becomes 
non-linear. Let’s consider, for example, BD gravity with non-linear interactions terms 

{3 + 2u}BD)\k‘^T =-^^Gpmd (150) 
in a Fonrier space. Here the interaction term X can be expanded as 

I{(p) = Mi{k)(p^ j - ki 2 )M 2 {ki,k 2 )(p{ki)(p{k 2 ) 

+ g J -- 5D{k-ki2fi)Mfiki,k2,kfipp{ki)ip{k2)(p{kfi), (151) 


where kij = ki + kj and kijk = ki + kj kk. 

The potential T is conpled to S through the BD scalar ip in a. fully non-linear way 
dne to the interaction term X. To derive the closed equations for 5m and 6 *^, we must 
employ the perturbative approach to Eq. (150). By solving Eq. (150) perturbatively 
assnming 99 < 1 , 4/ can be expressed in terms of 5 as [266] 




1 + X 


1 {k/a) 


3 U{k) _ 


Uk) + \{-] S(k), 


where 


k^ 


L1(A;) — X ( (3 + 2ci;bd)^ + Mi 


(152) 


(153) 


and = SttG. The fnnction S{k) is the non-linear source term which is obtained 
pertnrbatively as 


S{k) = 


Pmf 


d^kid^k2 , Jmiki) 5m{k2) 

5D{k — ki2)M2{ki, k2) 


54n(fc) J filnf 
(«;^Pm)^ f dz'kidz’k2dz’k^ 

486n(fc) J (27rp 

M2(fci, k2 + kfi}M2{k2, kfi) '1 5m{,ki) 5m{,k2)5m{,k^ 


n(fci)n(fc 2 ) 

5B{k - ki2fij\Mfiki,k2,kfii 


n(fc 


23 J 


n(fci)n(fc2)n(fc3 


( 154 ) 
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The expression (154) is valid up to the third-order in 6 m- Combining this with energy- 
momentum conservation, it is possible to calculate the density perturbation and the 
velocity divergence perturbatively and compute the non-linear corrections to the linear 
power spectrum (Fig. 10) [266, 267, 268]. 



Figure 10. The dark matter density-density, density-velocity, velocity-velocity power 
spectrum in ACDM (left) and f{R) gravity model with j/^ol = 10“^ (right). The 
solid lines are the results from the standard perturbation theory at the 1-loop level 
and symbols are results from simulations. From [267]. 


However, we cannot observe these power spectra directly. What we observe is the 
power spectrum of the galaxy over-density in the redshift space [269, 270]. There has 
been significant progress is developing non-linear models of the redshift space power 
spectrum (see for example [271] and references therein). In order to demonstrate the 
importance of modelling the non-linear effect in the redshift space power spectrum 
accurately, we consider the perturbation theory based template for the redshift power 
spectrum [272] 

P^^\k,fi)=DFoG[kR(TfiX |PKaiser(fc,fi) + ^(fc,/i) + 5(fc,h)}, (155) 

where /i = k^/k and the observer’s line-of-sight direction is taken to be the z-axis. The 
quantities pKaiser, and B are explicitly written as 


Tkaiser(fc,h) = P&s{k) - 2 Pse{k) + jo'^Peeik), 

A{k, /i) = —fc/i 


(156) 

k-p, -k) - B^{p, fc, -k-p)} ,(157) 


(27r)3 
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J M^F(p)F(k - p) ; (158) 

where Pss, Pee, and Pse respectively denote the anto-power spectra of the density and 
velocity divergence, and their cross power spectrnm. The fnnction is the cross 
bispectra defined by 

= (27r)^5£)(fci + ^2 + k^) k2, k^). (159) 

DpoG describes the damping of power dne to non-linear velocity dispersions known as 
the Finger of God effect 

Dpocikiaa^) = exp [-{kjaa^fi] . (160) 

We treat as a free parameter and marginalise over it when obtaining the constraint 
on model parameters. 

Fig. 11 demonstrates the importance of the non-linear corrections for an f{R) 
gravity model described by (24) [267]. We use the measured dark matter power spectrum 
in the redshift space at z = 1 from N-body simulations for f{R) gravity [273] as a mock 
observation. We now extract the model parameter \fRo\ by htting the redshift space 
power spectrum. We decompose the redshift power spectrum into the multipole power 
spectra 

Pf’(k) = J'dpP«>Uk, p)T,(p). (161) 

Fitting the measured monopole and the dipole P 2 ^ power spectra in simulations 
by the template (155) assuming that \fm\ and are free parameters, we study how well 
we can recover the input parameter |/i{o| in the simulation (in this case \fm\ = 10“^). If 
we do not include the A and B terms in (155), the recovered fjio is biased. If we include 
A and B but use the ACDM results, again the recovered fm is biased unless we cut off 
the power spectrum at low maximum wavenumber /Cmax < O.lh Mpc“^. Only the correct 
template reproduces the input parameter in an unbiased way up to /cmax = 0.15h Mpc“^. 
This implies that we need to specify the non-linear interactions in order to constrain the 
modihed gravity parameters accurately or it is required to use a very conservative cut-off 
in the wavenumber to remove the non-linear corrections. This is the major obstruction 
for model independent tests of gravity on linear scales. 

The same limitation applies to WL [274]. If we use a conservative cut-off in the 
analysis to remove the non-linear scales, the discriminatory power of weak lensing 
is signihcantly reduced. However, in order to include non-linear scales, we need to 
accurately model the non-linear power spectrum. For weak lensing, the perturbation 
theory is not adequate and N-body simulations are required. 
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Figure 11. Best-fit values of fuo as a function of the maximum wavenumber fcmax- We 
assume the cosmic variance limited survey of the volume V = 10/i“^Gpc^, we fit the 
PT template to the N-body simulation with |//jo| = 10“"* at z = 1. Filled circles are 
the results based on (155) in f{R) gravity, while filled triangles are the cases ignoring 
A and B terms. Open circles represents the results using A and B terms calculated in 
GR while crosses are the results ignoring the A and B terms and the damping term 
DfoG- From [267]. 


4 . 5 . Non-linear scales 

In order to find fully non-linear clustering of dark matter, it is required to perform N- 
body simulations even in the standard ACDM model. In GR, the Poission equation is 
linear so it is possible to superpose the forces. This was used to speed up the calculations 
by separating the long-range and short-range forces. In modihed gravity theories, the 
scalar held satishes the non-linear equation. Thus it is no longer possible to superpose 
the forces. The non-linear scalar held equation needs to be solved directly. 

There have been signihcant progress in developing N-body simulations [275, 276, 
277, 278, 279, 280, 281, 282, 283, 283, 284, 285, 286, 287, 288, 289, 290, 291, 292, 293, 
294, 295, 296, 297, 298]. Particularly, ECOSMOG [284], ISIS [291] and MG-GADET [296] solve 
a non-linear scalar held equation on a mesh in N-body simulations using the adaptive 
mesh technique (see Ref. [299] for a comparison between these codes). The ECOSMOG 
code for example has the following properties 

(i) It solves the scalar held on a mesh using the Newton-GauB-Seidal nonlinear 








Cosmological Tests of Modified Gravity 


47 


relaxation, which has good convergence properties. The density held on the mesh 
is obtained by assigning particles following the triangular-shaped-cloud scheme. 

(ii) The mesh can be adaptively rehned in high-density regions to achieve higher 
resolution and accuracy there, without affecting the overall performance. In order 
to speed up the convergence, a multigrid technique is used. 

(iii) It is efficiently parallelised, which makes the simulations fast. 



Figure 12. The relative difference of the power spectra of f{R) and ACDM models 
at z = 0. See the text for details. From [300] 

Fig.12 shows the power spectrum deviations from ACDM in the f{R) gravity model 
(24) with three different fm (F4: |/m| = 10"'^, F5: |/ro| = 10■^ F6 : |/ijo| = 10"® ) 
using various box sizes and resolutions obtained by ECOSMOG [300]. The chameleon 
mechanism works better for small |/ro| thus the power spectrum deviation from ACDM 
is more suppressed in F6. The dotted line is the linear theory prediction and the 
solid line is the prediction obtained by using the mapping formula {halofit [301]). The 
haloht model provides an accurate htting formula for the non-linear matter power 
spectrum obtained from a suite of N-body simulations in the ACDM model. It provides 
a prediction for the non-linear power spectrum for a given linear power spectrum. The 
haloht works better in F4 but as soon as the screening becomes important it fails to 
predict the the power spectrum deviation at high wavenumber (on small scales). This is 
not surprising as the haloht is calibrated in the ACDM model and it does not capture the 
ehect of the screening mechanism suppressing the deviation from ACDM on small scales. 
See Ref. [302] for a modihcation of halofit to account for the screening mechanism. This 
suppression of the power spectrum deviation from ACDM by the screening mechanism 
limits observational constraints. An interesting method has been proposed to enhance 
modihed gravity ehects by suppressing the contribution of the screened high-density 
regions in the matter power spectrum [303]. 

On small scales, the ehects from baryons become important. Hydrodynamical 
simulations have been developed for f{R) gravity and symmetron models Ref. [296, 
298, 293] but the degeneracy between modihed gravity ehects and baryonic physics 
remains to be understood. 
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N-body simulations enable us to understand how the screening mechanisms studied 
in section 3 work in the large scale structure of the universe. In the next section, we 
will look at the effect of the screening mechanisms on the structure formation in the 
Universe based on N-body simulations. 

5. Screening mechanisms in large scale structnre 

In this section, we consider two screening mechanisms, the chameleon and Vainshtein 
mechanism that we discussed in section 3 and study how GR is recovered in the 
large scale structure of the Universe and how we can distinguish between the different 
screening mechanisms. This section is based largely on Ref. [304]. 

The screening mechanisms are distinguished by how screened bodies fall in external 
helds [305]. As a consequence of universal coupling, all un-self-screened test bodies 
fall in the same way and obey a microscopic equivalence principle. In the chameleon 
and symmetron models, screened bodies do not respond to external helds while in the 
Vainshtein mechanism they do, as long as those helds have wavelengths long compared to 
the Vainshtein radius [306]. These diherences arise because of the non-superimposability 
of held solutions. Also, the peculiar structure of the non-linear interactions in the 
Vainshtein mechanism implies that the screening depends on the dimensionality of 
the system. For example, it does not work at all in one-dimensional systems. As a 
representative model for screening mechanism, we use the f{R) gravity model (24) that 
includes the chameleon model and the normal branch DGP model (27) that exhibits the 
Vainshtein mechanism. 

5.1. Screening of dark matter halos 

It is useful to construct analytic approximations for the scalar held solution inside dark 
matter halos. See Refs. [277, 307, 308, 309, 310, 311, 312, 313, 314, 315, 304] for analytic 
studies of dark matter halos in models with the above screening mechanisms. Here we 
mainly follow Refs. [308, 315, 304]. 

We assume the dark matter density prohle is described by the NFW prohle [316] 
with a mass M^. The mass is dehned as the mass contained within the radius r = r^,; 
the density at ta is PcutA, where pcut is the critical density of the Universe. The NFW 
prohle is given by 

p{r) = p.f ()-) . m = (162) 

where p* = p{rs) is hxed so that the mass within ta is Ma. The scale radius is 
more conveniently parameterised by the concentration ca = r^/rg. By integrating this 
density prohle, we obtain the enclosed mass within the radius r, M(< r), as 

= F(y) = -fi—+h-,(l + y). (163) 

F[c^) l + y 
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The scalar field equation in the nDGP model, Eq. (30), can be solved analytically 


[308] 

dip _ GM{< r) 4 
dr 3/3' 

where r* is the Vainshtein radius 


^(^) ’ 9{x) = , 


r* = 


/ 16GM(< r)r 


I 


9/32 


1/3 


(164) 


(165) 


Inside the Vainshtein radius, the scalar force is suppressed compared with the Newtonian 
potential due to the non-linear derivative interactions. Outside the Vainshtein radius, 
the linear solution is realised where g{r/rfi) —)■ 1/2. For a larger Tc, the Vainshtein radius 
is larger thus the region in which the fifth force is suppressed becomes larger and we 
recover GR. 

The scalar field in f{R) gravity can be approximated as [308] 

lG(M(<r)-M(<r.„)) 

dr 


(166) 


3 r2 

for r > Tscr and dp/dr = 0 for r < r^cr if the thin shell condition (100) is satisfied. The 
screening radius for the NEW profile is obtained as [317] 


_ /2 UnKj-a) : 
' 3 |/«„|F(ci) 


(167) 


where the Newtonian potential is given by |4/7v| = GM^/r^. The screening radius is 
determined by the ratio between the Newtonian potential of the halo and \fm\, thus it 
is mass dependent. 

We define the ratio between the fifth force and the Newton force as 

dT jv 


Am - jr, 


F - 

® 2 dr’ 


Fg = 


dr 


(168) 


For linear solutions without screening. Am = 1/3/3 in nDGP and Am = 1/3 in f{R). 


5.2. Simulations 

We consider iV-body simulations for nDGP and f{R) models performed by the 
EGOSMOG code described in section 4.5 [304]. To highlight the difference in the 
screening mechanisms, we simulate pairs of nDGP and f{R) models with the AGDM 
background and identical a.t z = 0, which roughly removes the difference in screening 
due to the difference in linear growth. Specifically, we simulate three f{R) models: F4 
(I/foI = 10“^), F5 (I/hoI = 10"®), and F6 (|/ro| = 10"®), and three corresponding 
nDGP models whose parameters are listed in Table I. 

The simulations are visualised in Fig 13, where the projected density field for AGDM 
and two modified gravity models are shown at z = 0 [304]. As shown, the structures are 
more clustered in nDGP and f{R) models due to the enhanced gravity, although the 
enhancement of the clustering in these models are different. To quantify the difference in 
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fiR) 

nDGP 


F4: \fno\ = lO"" 

nDGPl: HoTc = 0.57 

as = 0.946 

F5: \fno\ = 10"^ 

nDGP2: ifoG = 1-20 

as = 0.891 

F6: \fm\ = IQ-® 

nDGP3: HoVc = 5.65 

as = 0.854 


Table 1. The parameters for f{R) and nDGP simulations. The baseline cosmology 
was chosen to be the model favoured by recent Planck observations [18]: = 

0.022161, = 0.11889, = 0, h = 0.6777, = 0.9611, and erg = 0.841. The 

simulations use 256^ particles in a L = 64/i“^ Mpe box from the initial redshift z = 49 
down to z = 0. Each set of models {f{R), nDGP, and AGDM) are simulated using the 
same initial condition, which is generated using the Zeldovich approximation, and we 
simulate three realisations for each model to reduce the sample variance. See [304] for 
details. 


LCDM nDGP F4 



y (Mpe) y (Mpe) y (Mpe) 


Figure 13. The snapshots (64Mpe x 64 h~^ Mpe) showing the projected density 
field at z = 0 for LCDM (left), nDGPl (middle) and the F4 (right) models. From [304] 
(published on 29 July 2015 ©SISSA Medialab Sri. Reproduced by permission of lOP 
Publishing. All rights reserved). 


clustering, we show the fractional difference of the power spectrum in f{R) and nDGP 
models with respect to ACDM at z = 0 in Fig. 14 [304]. The dotted lines show the 
linear prediction while the dash-dotted lines show the halofit prediction [301]. 

The agreement of the linear prediction with simulations is better for the nDGP 
model. The excellent agreement with linear theory on large scales is one of the key 
features of the Vainshtein mechanism. This is because when the Vainshtein screening 
mechanism works, even if the hfth force is suppressed inside halos, these ‘screened’ halos 
can still feel external scalar helds as long as those helds have wavelengths longer than 
the Vainshtein radius [306]. We will conhrm this picture later by studying the velocities 
of dark matter particles inside halos. On small scales, the Vainshtein mechanism is very 
effective and the power spectrum deviation approaches zero quickly at high k. 

On the other hand, in the f{R) model with the chameleon screening mechanism, 
once the hfth force is suppressed inside dark matter halos, these screened halos no 
longer feel the external hfth force invalidating the use of the linear theory even on large 
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k [h/Mpc] 

Figure 14. The fractional difference in power spectrum for nDGP (left) and f{R) 
(right) models with respect to the ACDM model. The data points with error bars 
show the simulation result, and the dashed (dash-dotted) curves show the Haloflt 
(linear) predictions. From [304] (published on 29 July 2015 ©SISSA Medialab Sri. 
Reproduced by permission of lOP Publishing. All rights reserved). 


scales. The chameleon mechanism is not as efficient as the Vainshtein mechanism. Thus 
the power spectrum deviation does not approach zero even though it is signihcantly 
suppressed compared with the linear prediction. 


5.3. Screening of dark matter particles 


We hrst study how screening operates at the level of particles [315, 304]. The cosmic 
web of large scale structure consists of an interconnected hierarchy of halos, filaments, 
walls, and voids. Here we use ORIGAMI to determine the cosmic web morphology 
of each dark matter particle in a simulation, which compares final positions to initial 
Lagrangian positions to determine whether a particle has undergone shell-crossing along 
a given set of axes (see Ref. [318] for details). Shell-crossing denotes the formation of 
caustics within which the velocity field is multi-valued, called the multi-stream regime. 
The number of orthogonal axes along which shell-crossing has occurred corresponds to 
the particle’s cosmic web morphology and denoted by the morphology index M: halo 
particles have crossed along three axes (M = 3), hlaments along two (M = 2), walls 
along one (M = 1), and void particles are in the single-stream regime (M = 0). 

We quantify the screening by calculating the deviation of the fifth force to 
gravitational force ratio from the linear relation {Am = 1/3/3 in nDGP and Am = 1/3 
in f{K)) 


AF 


^5 

AmFg 


- 1 , 


(169) 


which ranges between AF = — 1 when screening is working to AF = 0 when it is not. 
For some particles AF can be greater than 0 due to numerical noise, especially for low 
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values of F 5 and Fq (see, e.g., [319]). 

In Figure 15 we show histograms of AF split according to ORIGAMI morphology 
for each modihed gravity model [304], Note that the histograms are normalised to 
peak at unity so that the shape of each can be seen, so they do not reflect the relative 
abundances of the particles according to their cosmic web morphology. It is clear that, as 
for all 3 nDGP models the halo particles are screened while the hlament, void, and wall 
particles are unscreened. This reflects the dimensionality dependence of the Vainshtein 
mechanism [181, 320, 315]. The non-linear term in Eq. (30), — {didjipfi, vanishes 

for a one-dimensional system. The distribution of AF for halo particles is still broad, 
reflecting the fact that screening becomes weaker beyond the virial radius, so particles 
in the halos’ outer edges can have larger values of AF. 



Figure 15. Cosmic web morphology dependence of the Vainshtein (left) and 
chameleon (right) screening mechanisms, given by histograms of AF, the deviation 
from the linear relation of the fifth force to gravitational force ratio. In DGP models 
with the Vainshtein mechanism, there is a clear difference between halo particles (solid, 
red line) and filament, wall, and void particles, while no such distinction exists for f{R) 
models with the chameleon screening mechanism. From [304] (published on 29 July 
2015 ©SISSA Medialab Sri. Reproduced by permission of lOP Publishing. All rights 
reserved). 

On the other hand, there is no such morphology dependence of the screening 
mechanism for the chameleon models on the right side of Figure 15. For the halos. 
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screening is not very effective in F4 and most particles follow the linear relation; in F5 the 
halo particles have a small peak at AF = — 1 as the most massive halos become screened; 
and in F6 most of the halo particles are in screened halos, while some remain unscreened. 
Note that since the histograms are for halo particles and not halos themselves, massive 
halos are weighted more heavily, resulting in many screened halo particles in F6, while 
the smaller bump of unscreened halo particles is due to halos with low mass. We will look 
at the screening of halos in the next section. The wall and void histograms notably do 
not peak at AF = 0 in F4, and the hfth force is further suppressed in F5 and especially 
F6. This is because the background Compton wavelength is quite short, ~ IMpc in F6, 
and the scalar held does not propagate beyond this length, providing a blanket screening 
for particles that are sparsely distributed. The hlament distribution develops a double 
peak in F6: a narrow peak of unscreened hlament particles, which have large forces and 
are in relatively dense environments, and a broader peak of low AF hlament particles 
that are blanket screened. 

5 . 4 . Dark Matter Halos 

We now turn to a comparison of the screening of halos in the Vainshtein and chameleon 
mechanisms [304]. 

5 . 4 . 1 . Screening - mass dependence As with the dark matter particles, to determine 
whether screening is ehective we calculate the ratio of the hfth force to gravitational 
force. Am- The value of Am for each halo is given by the average Am of all the particles 
in the halo within the halo’s virial radius, i? 2 oo- This is plotted as a function of halo 
mass, M 200 , for all nDGP and f{R) simulations in Figure 16 [304]. 
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Figure 16. Ratio of the fifth force to gravitational force as a function of halo mass for 
nDGP (left) and f{R) (right) models. From [304] (published on 29 July 2015 ©SISSA 
Medialab Sri. Reproduced by permission of lOP Publishing. All rights reserved). 

For the DGP models, the Am of the halos depends on the model parameter and is 
independent of mass [308, 315]; as the model parameter changes to make the deviation 
from ACDM stronger. Am increases. However, note that all halos are screened in the 
Vainshtein mechanism: the linear values of Am for nDGPl, nDGP2, and nDGP3 are 
















Cosmological Tests of Modified Gravity 


54 


0.20, 0.11, and 0.03, respectively, well above the corresponding values in Figure 16. As 
we will see in the next section, the Vainshtein suppression gradually reduces (and Am 
increases) outside the virial radius, but even including these particles in the calculation 
of Am only increases Am by ~ 50 — 70% and the halos remain screened overall. 

In contrast to the Vainshtein mechanism. Figure 16 shows there is a clear 
dependence on both mass and model parameter in the chameleon mechanism [308, 304]. 
When the deviation from ACDM is high, in F4, screening becomes ineffective for all 
halos; in F5, screening is effective only for high mass halos; and in F6, there is a 
population of unscreened small halos and a transition to screened large halos. Including 
particles outside the virial radius in the calculation of Am has a very small effect, 
increasing Am for some halos but not changing the overall trends; we will show in the 
next section that the radius of transition from screened to unscreened parts of the halo 
depends on the halo mass in the chameleon mechanism. 


nDGP2 


rDGP3 



Figure 17. Profiles of Am in logarithmic bins of normalised radius, for nDGP2 (left 
panel) and nDGP3 (right panel), with the analytic prediction for an NFW profile 
plotted in red. Screening profiles are independent of halo mass, and Vainshtein 
screening suppresses the fifth force within the virial radius. From [304] (published 
on 29 July 2015 ©SISSA Medialab Sri. Reproduced by permission of lOP Publishing. 
All rights reserved). 

5 . 4 . 2 . Screening profiles To determine the radial dependence of the screening within 
the halos, here we calculate the ratio between the fifth force and Newtonian force 
as a function of normalised radius R/R200 by averaging Am of the halo particles in 
logarithmic radial bins. Am profiles are plotted in Figure 17 for nDGP2 and nDGP3 
[304]. It is clear that regardless of halo mass, the Vainshtein screening profiles of dark 
matter halos are roughly the same and correspond quite well to the spherically symmetric 
analytic solution for an NFW profile (164) [315]. For both models, the fifth force is 
mostly suppressed within the virial radius and Am increases outside the virial radius, 
and the magnitude of this increase is greater for models with a stronger enhancement 
to gravity, nDPG2 (and nDGPl, not shown). 

In f{R), Figure 16 shows that chameleon screening does depend on mass, so in 
Figure 18 we split up the Am profiles into four different mass bins for the F6 model 
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Figure 18. Profiles of Am in logarithmic bins of normalised radius for the F6 model, 
split into four bins of halo mass, with the analytic prediction calculated for each mass 
bin plotted in red. Screening prohles in the chameleon mechanism depend on halo 
mass. From [304] (published on 29 July 2015 ©SISSA Medialab Sri. Reproduced by 
permission of lOP Publishing. All rights reserved). 


[304]. There is a much sharper and less gradual transition from the chameleon screened 
inner regions of the halo to the unscreened outer regions compared to the Vainshtein 
mechanism, and the radius of this transition depends on the halo mass. There is again 
quite a good agreement with the analytic predictions Eq. (166) [308, 304]. 

Both Figure 17 and Figure 18 highlight the importance of probing galactic halos in 
unscreened regions beyond their virial radii in order to detect deviations from ACDM, 
for all halo masses if the Vainshtein mechanism is operating and for high mass halos in 
the chameleon mechanism. For example. Refs .[321, 322] proposed a method to measure 
the velocity held by stacking redshifts of surround galaxies around galaxy clusters from 
a spectroscopic sample. It was shown that order unity deviations from ACDM show 
up both in f{R) gravity and the nDGP model. On the other hand, the effect on the 
interior mass prohle, which is measurable through stacked weak lensing is much less 
affected by a modihcation of gravity (see also [323]). Ref. [324] proposed to use the 
2D galaxy velocity distribution in the cluster infall region by applying the galaxy infall 
kinetics model developed by Ref. [325]. 

5. 5. Velocities 

The non-linear derivative interaction that is responsible for the Vainshtein mechanism 
enjoys the Galilean symmetry (33), 9^0 —)■ Cfj,, which means that if the external 
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fields have wavelengths that are long compared to the Vainshtein radius, Eq. (165), we 
can regard the gradient of these external helds as constant gradients in the vicinity of 
an object. We can always add these constant gradients to the internal held generated by 
the object, and thus the internal held can superimpose with external helds. Even if the 
internal held is suppressed by the Vainshtein mechanism, the object still feels the hfth 
force generated by the external helds [306]. On the other hand, in screening mechanisms 
that rely on non-linearity in the potential or the coupling function to matter, such as 
the chameleon and symmetron mechanisms, the internal held generated by an object 
does not superimpose with an external held. Therefore, the held inside the object loses 
knowledge of any exterior gradient and the hfth force generated by the external held, 
and thus once the object is screened, it does not feel any hfth force [305, 167, 306]. 

Velocity Dispersion Peculior Velocity 




Figure 19. Velocity dispersion (left) and peculiar velocity (right) ratios of nDGP (top) 
and f{R) (bottom) models with respect to matched halos in the ACDM simulations, as 
a function of halo mass. Lines indicate the average value in each mass bin, and error 
bars indicate the scatter. From [304] (published on 29 July 2015 ©SISSA Medialab 
Sri. Reproduced by permission of lOP Publishing. All rights reserved). 

In Figure 19, we plot the velocity dispersion ratios and peculiar velocity ratios for 
all DGP and f{R) models [304]. The velocity dispersion ratios are scaled by the virial 
expectation, oc to remove the standard mass dependence [308]. The lines are 

the average values in bins of mass, and the error bars show the Icr standard deviation. 
Both the velocity dispersion and peculiar velocity ratios show no dependence on mass 
for the Vainshtein mechanism [308, 315], but the peculiar velocity ratios deviate from 
0, especially in the nDGPl model, due to the effect of external helds. 
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The f{R) models, on the other hand, do show a mass dependence in the velocity 
dispersion ratios: in F4, where screening is not effective, the velocity dispersion is 
enhanced; in F5, the velocity dispersion is only suppressed at high masses; and in F6, 
the transition between enhanced and suppressed velocity dispersion occurs at an even 
lower mass. These trends mimic the mass dependence of the screening ratio. Am, in 
Figure 16. However, the peculiar velocity ratios do not show this mass dependence; they 
are suppressed in F6, somewhat enhanced in F5, and further enhanced in F4. Since most 
of the particles in a halo are located near the centre, the peculiar velocity effectively 
probes the halo centre and is suppressed because the halo core is screened, as seen in 
the Am profiles for F6 in Figure 18. Unlike in the Vainshtein mechanism, once screened 
by the chameleon mechanism a halo does not feel the effect of external fields and so its 
peculiar velocity is suppressed. 

This casts doubt on the effectiveness of using redshift space distortions to detect 
f{R) gravity at the level of F6. The Vainshtein mechanism is better suited to tests in 
the linear regime, since even though halos are screened, they can still feel the effect of 
external fields induced by large scale structure. 

5.6. Screening - environment dependence 

In the previous sections, we have shown that chameleon screening depends on both the 
mass of the object and the model parameter, while Vainshtein screening is independent 
of halo mass [308]. In Section 5.3 we showed that Vainshtein screening of dark matter 
particles depends on their cosmic web morphology [315]; halo particles are screened while 
filament, void, and wall particles are not, reflecting the dimensionality dependence of the 
Vainshtein mechanism [181, 245]. However, it was found that the Vainshtein screening 
of halos themselves does not depend on their large scale cosmic web environment [315], 
and the chameleon mechanism has no cosmic web dependence for either particles or 
halos, so here we use a different definition of halo environment. 

We use a density-based definition of environment developed in Ref. [326], 

„ _ dN,MN/M>f 

L>NJ = -• 

rN 

This is the distance dj^ to the iVth nearest neighbour having mass, M^r, at least / 
times as large as the halo mass, M, scaled by the virial radius of the neighbouring 
halo, r^r. Dl l (hereafter just D) is almost uncorrelated with the halo mass [326], and 
several studies have found that both chameleon and symmetron screening mechanisms 
correlate with this environment parameter [327, 328, 329, 330, 319]. In particular, 
while massive halos can be self-screened and in general smaller halos can be unscreened, 
small halos that live in dense environments (where D is small) can be environmentally 
screened. Subhalos, especially those within the virial radius of their host halo with 
D < 1, are usually environmentally screened. We use the AHF halo finder [331] to 
study the environmental dependence of chameleon and Vainshtein screening [304]. 

In Figure 20 we show the results for chameleon screening [304]. In the left panel 
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we plot Am vs. halo mass, M^ir, with a logarithmic scaling of Am instead of the 
linear scaling of ORIGAMI valnes in Fignre 16. The same trend in mass dependence 
of screening is seen, except AHF halos can have mnch lower valnes of Am, which 
is dne to the presence of snbhalos. In the right panel of Figure 20 we show the 
environmental dependence of chameleon screening, after first removing halos with mass 
above 10^^ h~^ Mq that are self-screened. Halos in dense environments, and especially 
those with D < 1, tend to be screened while those in under dense environments, 
especially those with D > 10, are unscreened [327]. 



Figure 20. The ratio of the hfth force to gravitational force for AHF halos in f{R) 
models as a function of mass (left panel) and environment D (right panel). Only halos 
with mass below 10^^ h~^ Mq are shown in the right panel to remove halos that are 
massive enough to be self-screened. Chameleon screening displays both a mass and 
environmental dependence. From [304] (published on 29 July 2015 ©SISSA Medialab 
Sri. Reproduced by permission of lOP Publishing. All rights reserved). 

We show the profile of the mass difference as a function of rescaled radius for the 
l/itol = model in Fig 21 [327]. To see the environmental effect on the prohle, we 
show the result for the samples selected according to both their mass and D. As we can 
see, the small halos in the under dense region are almost not screened at any radius, 
while the halos with similar mass in the dense region are efficiently screened, and the 
screening effect is stronger in the core of the halos. For the large halos, the innermost 
part is well screened regardless of external environment due to the high matter density 
there, but the part close to the edge shows a clear environmental dependence, and the 
difference can be as large as 3 orders of magnitude in Am in different environments. 
This is because in this region the external environment plays an important role. 

The results for the Vainshtein screening mechanism are shown in Figure 22 [304]. 
The left panel shows the screening as a function of mass, with again no mass dependence, 
though there is a lot of scatter in Am for low mass halos. The right panel shows Am 
as a function of environment: in general there is no environmental dependence, but 
halos within their host halo with D < 1 are screened even further by their host halo. 
Note again that all halos are screened in the Vainshtein mechanism: the linear values 
oi Am for nDGPl, nDGP2, and nDGP3 are 0.20, 0.11, and 0.03, respectively, well 
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Figure 21. The profile of logioAM as a function of the rescaled halo radius r/r 34 o 
for the |/flo| = 10“® model. We show the profile with 1 — a error bars for the halos 
divided into four categories as illustrated in the legend. The red dashed line shows 
Am(?’34o) = 1/3. From [327]. 



Figure 22. The ratio of the fifth force to gravitational force for AHF halos in nDGP 
models as a function of mass (left panel) and environment D (right panel). There is 
no mass or environmental dependence of screening in the Vainshtein mechanism, and 
sub-halos within their host virial radius (with D < 1) have fifth forces that are even 
further suppressed. From [304] (published on 29 July 2015 ©SISSA Medialab Sri. 
Reproduced by permission of lOP Publishing. All rights reserved). 


above the corresponding values in Figure 22. The Vainshtein mechanism is therefore 
very efficient at screening halos: there is no dependence on mass [308, 315], cosmic web 
morphology [315], or density of their local environment. 
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The character of the screening mechanism leaves distinct signatures in non-linear 
structures as we discussed in the previous section. In this section we review novel 
astrophysical tests of gravity that have been developed recently. These tests are model 
dependent as opposed to the model independent tests of gravity on linear scales. We 
need to specify the screening mechanism to discuss the tests. 

6.1. Chameleon mechanisms 

We hrst consider the chameleon mechanism. The chameleon screening of dark matter 
halos depends on their mass and environment. Thus we need to carefully hnd the places 
where gravity is modihed. From Fig. 21, the places where gravity is modihed are at 
the outskirt of large dark matter halos and small dark matter halos in the under dense 
environment. 

6.1.1. Clusters The abundance of clusters already gives strong constraints on |/ijo| 
[332, 333, 334]. For large Ifml > 10“® the screening is ineffective and stronger gravity 
creates much more massive dark matter halos as is seen in Fig. 23. On the other hand, 
for Ifml < 10“^ the screening is effective for massive halos and the relative enhancement 
of the mass function disappears. The latest constraints on |/^o| from cluster abundance 
is obtained as Ifml < 2.6 x 10“® for n = 1 in Eq. (24) [334]. 

If the screening length is smaller than the radius of the dark matter halo, we will see 
the transition from the screened and unscreened region through the difference between 
the lensing and dynamical masses [317] . We assume that a cluster of galaxies satisfies 
hydrostatic equilibrium 

dPtehrmal _ Giff(< r) fi d(f) flTlI 

dr r Mpi dr ’ 

where Pthermai is the thermal pressure of the gas. The enclosed mass M(< r) agrees 
with the lensing mass as the scalar field does not affect the lensing. The last term in 
Eq. (171) is the contribution from the fifth force. This changes the dynamical mass, 
which can be inferred from the pressure gradient of the gas. Fig. 24 shows the enclosed 
mass reconstructed from X-ray observations using the hydrostatic equilibrium equation 
with the chameleon force for the Coma cluster [335]. It also shows the measurement 
of the lensing mass for the Coma cluster [336]. Since the chameleon gravity adds an 
additional contribution to the dynamical mass, the lensing mass estimated from the X- 
ray observation decreases. In this way, we can put constraints on the chameleon gravity. 
Using the Coma cluster [335] and the stacked 58 clusters between 0.1 < x < 1.2 using 
X-ray data from XMM Cluster Survey and weak lensing data from the CFHTLenS [337], 
the constraint on f{R) gravity was obtained Ifml < 6 x 10“®. 

6.1.2. Dwarf galaxies As is clear from Fig. 20 and Fig. 23, for |/^o| < 10“®, cluster are 
effectively screened and we need to find objects with a shallower potential. Another place 
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Figure 23. The relative enhancement of the halo mass function in f{R) gravity 
models compared with ACDM. The solid line and dashed line are analytic predictions 
(see [312] for the detail). From [312]. 


where gravity is strongly modified is provided by dwarf galaxies in voids. Dwarf galaxies 
have a potential ~ 10“^ thus they are not self-screened even with Ifaol = 10“®. 
We need to carefully identity dwarf galaxies that are not environmentally screened. 
One way is to use the D parameter introduced in the previous section (Eq. (170)). 
Another environmental measure is defined in terms of the sum of gravitation potentials 
of neighbouring galaxies [328] 


$ext 


di<\c+ri 


GMj 

di 


( 172 ) 


where di is the 3D distance to the neighbouring galaxy with mass and the virial 
radius r*, Ac is the background Compton wavelength of the scalar field. Using the thin 
shell condition (100), we select the environmentally unscreened galaxies 

"hext < 2 I/™I- (I’^S) 

This parameter is more suited to identify the environment of galaxies when uncertainties 
in the estimation of mass and distance and the incompleteness of the samples are 
included compared with D [328]. Using this criteria a screening map has been created 
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Figure 24. Radial mass profile of the Coma cluster. The shaded region is the 
observationally allowed Icr region of the WL observation. The blue solid curve is the 
thermal mass component estimated from the X-ray observations. The red solid and 
dashed curves are the combination of the termal and chameleon mass components. 
From [335] (published on 14 April 2014 ©SISSA Medialab Sri. Reproduced by 
permission of lOP Publishing. All rights reserved). 


using SDSS galaxies within 200 h ^Mpc [328] (see Fig. 25 and 26). We also introduce 
the self-unscreening condition 

‘f’int < (174) 

where <hint = GM/R, M is the mass and R is the virial radius of a halo. 

There are a number of effects that we expect to see in unscreened galaxies 
summarised in [338]. 

• Offset between stellar and gaseous component. 

In the presence of the external scalar held, unscreened components such as dark 
matter or neutral hydrogen (HI) will fall faster along the external held than will 
stars that are self-screened. This diherence may cause an ohset between the centroid 
of the stellar disk and the gaseous disk. 

• Warping of the stellar disk. 

As the dark matter halo moves along the external potential gradient, the screened 
stellar component lags behind the dark matter halo. If the potential gradient is 
aligned with the axis of rotation, this may introduce a U-shaped warp in the stellar 
disk. The warp is expected to align with the potential gradient. 

• Asymmetry in stellar rotation curve. 

The ohsets between stellar and halos centres can perturb the stellar disk and cause 
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Figure 25. Scatter plot for <i>ext vs. the dynamical mass Mdyn (equivalent to |$int|), 
coloured accordingly to the difference between the dynamical and lensing mass Am- 
A simple cut in $ext and Mdyn given by Eqs. (173) and (174) seem to work to separate 
screened from unscreened galaxies. From [328] (18 July 2012 ©SISSA Medialab Sri. 
Reproduced by permission of lOP Publishing. All rights reserved). 



Figure 26. The environmental screening map generated in the SDSS region at 210Mpc 
for 10 X 10 degree of survey area, approximately 38 Mpc, with 2 arcmin resolution. The 
external potential <i>ext is shown with the screening condition evaluated using Eq. (173) 
for models with |/i^o| = 10“^,10“® and 10“^ (left, middle and right panels). The cut 
in screening classification 3/2|/ijo| is shown in the map as a white contour line (regions 
inside the contour are screened). From [328] (18 July 2012 ©SISSA Medialab Sri. 
Reproduced by permission of lOP Publishing. All rights reserved). 


an asymmetry in the stellar rotation curve as the dominant force on the stellar disk 
is from the potential of dark matter halos. 

• Mismatch of the rotation velocities of a gaseous and stellar disk 

The rotation velocity of an unscreened gaseous disk is enhanced relative to the 
screened stellar disk by a factor of ^/ImKgJG. This test requires us to derive the 
stellar rotation curves from stellar absorption lines instead of using Ha or the 21 
cm line as these lines prove the unscreened gaseous components of the galaxies. 
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All these tests have been carried out recently [338, 339]. The absence of these effects 
give the strong constraint |/ijo| < 10“®. These constraints can be improved by on-going 
surveys such as Sloan Digital Sky Survey (SDSS)IV Mapping Nearby Galaxies at APO 
(MaNGA) which will hnd many unscreened galaxies even if Ifaol = 10“® [340]. 

There are also interesting implications for dwarf galaxies in the Milky way if the 
screening radius rscr of the Milky way is 60 kpc [341]. 

6.1.3. Stars Another astrophysical object which has a shallow potential is a post main- 
sequence star whose potential is |T 7 v| ~ 10“^. The stellar structure is again determined 
by the hydrostatic equilibrium. The effect of the hfth force can be approximated by 
Eq. (166). Due to the enhanced gravity, stars are expected be more luminous [342, 343]. 
This modihed equation has been implemented in the stellar structure code MESA [344]. 
Assuming that a star is in an unscreened galaxy, the evolution of a 1 solar mass star 
in the Hertzprung-Russell diagram is shown in Fig. 27 in the case of f{R) gravity. For 
Ifml = 0.67 X 10“^, the main-sequence stars are screened and the modihed gravity 
affects only red giants. For |/ijo| = 0.67 x 10“® and 3.35 x 10“®, main sequence stars 
are also affected. 



logTgff 


Figure 27. The Hertzprung-Russell diagram for stars of one solar mass with initial 
metallicity Z = 0.02 in f{R) gravity. The black line shows the tracks for stars in GR 
while the red, bine and green tracks correspond to stars in modified gravity models 
with the chameleon mechanism with \ fRo\ = 0.67 x 10“^, 0.67 x 10“® and 3.35 x 10“®. 
From [343]. 

This modihcation of the stellar structure was used to put a strong constraint on 
the chameleon mechanism [345]. Gepheid variable stars with 5 — IOMq pulsate with a 
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known period-luminosity relation once they go off the main sequence and situated in a 
narrow temperature gap in the HR diagram. This property has been used to measure 
the distance to the Cepheids. The period r is determined by the strength of gravity 
r oc . This relation is calibrated using local group stars that are screened. Hence 
if one uses this relation to estimate the distance, then the derived distance will be 
incorrect. If gravity is stronger than GR, the period of Cepheids becomes shorter. It is 
found that the change in the distance is given by 


M 

If 


-0.3 


AG 

~G~' 


( 175 ) 


If we measure the distance using Cepheid assuming GR, the measured distance will 
be shorter than the actual distance. In order to test this, we need another distance 
indicator that is insensitive to gravity and measure the correct distance. The tip of the 
red giant branch (TRGB) distance is such a distance indicator as the absolute magnitude 
of red giants that go off the red giant branch after a helium flash is insensitive to 
gravity. The comparison between the two distances has been performed and compared 
against the prediction of the chameleon mechanisms. This led to the stringent constraint 
I/roI < 5 X 10-h 


6.1.4- Summary of constraints Fig. 28 summarises the constraints on chameleon 
gravity in the case of f{R) gravity. A detailed summary of constraints is available 
in Ref [346]. 

Gonstraints from cosmological probes such as GMB and large scale structure are 
relatively weak |/ijo| < 10“^ — 10“^. Glusters give an order of magnitude stronger 
constraint Ifml < 10“^ which requires the understanding of non-linear structure 
formation. In order to satisfy the Solar System constraints, the Milky Way needs to be 
screened. This condition depends on whether the Local Group gives the environmental 
screening and the constraint is given by |/ijo| < 10“"^ — 10“®. Gonstraints from dwarf 
galaxies and Gepheid are much stronger |/ijo| < 10“^. All the constraints except for the 
Solar System constraint have been obtained in the last 10 years illustrating impressive 
development of new tests of gravity on cosmological and astrophysical scales. 


6.2. Vainshtein mechanism 

Unlike the chameleon mechanism, it is very difficult to test models with the Vainshtein 
mechanism on small scales. The screening of dark matter halos does not depend on 
mass nor environment and all dark matter halos are screened inside the virial radius. In 
order to hnd modihcations of gravity, we need to go beyond the virial radius of halos. 

An interesting test using black holes at the centre of galaxies has been proposed 
[306]. As we have seen in section 5.5, even if a dark matter halo is “screened” and the 
hfth force is suppressed within the dark matter halo, it still feels the external gradient 
of the scalar held as it still carries a scalar charge. On the other hand, black holes do 
not carry a scalar charge due to the no-hair theorem. Thus the central black hole will 
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Figure 28. Astrophysical and cosmological limits on chameleon theories. The spatial 
scale on the x-axis gives the range of the length scales probed by particular experiments. 
The parameter on the y-axis is the background field |//jo| value or the range of the 
interaction for an /(i?) gravity model. The rectangular region gives the exclusion zone. 
From [13]. 


lag behind the dark matter halos and there may be a displacement of the central black 
holes. The magnitude of this effect depends on the external held and the strength of 
the hfth force and it is estimated as [306] 


r = 0.1 kpc X 2/3^ 


IVT 


ext 


O.OlMpc 

Po 


-3 


( 176 ) 


20 (/cm/s)2/kpc^ 
where po is the central density of the galaxies. 

In beyond Horndeski theories discussed in section 2.5, the Vainshtein mechanism is 
broken inside matter. The Poisson equation within the Vainshtein radius is given by 


_0 _ GM{< r) T dfiM{< r) 

dr 4 

where T characterises the deviation from the Horndeski theory. The last term modihes 
GR even within the Vainshtein radius. This for example changes the structure of stars, 
rotation velocities of galaxies and lensing [347, 348]. 

Observational tests have not been carried out yet and there remains to be seen 
if we can get tighter constraints on the model parameter such as Uc compared with 
cosmological observations and the Solar System constraint from these astrophysical 
tests. 
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The discovery of the accelerated expansion of the Universe has come relatively late in onr 
stndy of the cosmos, bnt in showing that gravity can act repnlsively, it has opened up 
many new questions about the nature of gravity and what the Universe might contain. 
Is the acceleration being driven by dark energy? Or is general relativity itself in error, 
requiring a modihcation at large scales to account for the late acceleration? 

Modihed gravity models have provided many interesting new ideas to solve the 
cosmological constant problem by modifying the way in which the vacuum energy 
gravitates. The expansion of the Universe can accelerate without the cosmological 
constant due to the modihcation of gravity on cosmological scales. Gravitons may have 
a mass and provide an effective cosmological constant realising the technical naturalness 
of the smallness of the cosmological constant. Einstein equations are the only second- 
order local equations of motion for metric derivable from the action in 4D. If we modify 
GR, a new degree of freedom appears in the graviton, modifying gravity even in the Solar 
System. One of the most important recent developments are the screening mechanisms 
that enable us to modify gravity signihcantly on cosmological scales while satisfying the 
stringent Solar System constraints. Despite extensive study of modihed gravity models, 
we still do not have a model that can be an alternative to the AGDM model. Our quest 
to hnd modihed gravity models still continue. 

Structure formation in our Universe can be diherent even if the geometry of 
the homogeneous and isotropic universe is the same in these two classes of models, 
ohering a possibility to distinguish between them observationally. It is in principle 
possible to construct non-parametric consistency tests of GR on cosmological scales by 
combining various probes of large-scale structure, because the Einstein equations enforce 
a particular relation between observables. On large scales, there are two functions 
of time and space that are required to describe these relations. One describes the 
relation between the density perturbation and the Newton potential, which determines 
the motion of non-relativistic objects such as galaxies; the other parametrises the 
relation between the density perturbation and the lensing potential, which determines 
the geodesics of photons. In the next hve years, new surveys such as DES and eBOSS will 
be able to constrain several parameters at the 5-10% level. This shows that upcoming 
surveys will make cosmological tests of gravity a reality in the next hve years. Surveys 
such as Euclid will further make cosmological tests of gravity a precision test providing 
constraints on these parameters at the 1% level. 

It is possible to construct model independent tests on linear scales, but the 
bulk of information is available on non-linear scales. Also in order to extract the 
information of the linear growth rate of structure, it is required to model the non¬ 
linear corrections accurately. The non-linear structure formation is complicated by the 
screening mechanisms that restore GR on small scales. Perturbation theory and N-body 
simulations have enabled us to understand how structure formation is affected by the 
screening mechanisms. The screening mechanisms leave distinct signatures in the non- 
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linear structure formation. This has led to the developments of novel astrophysical tests 
of gravity. In the case of the chameleon mechanism, the astrophysical tests using dwarf 
galaxies and stars provide stringent constraints on the model that are even stronger 
than the Solar System constraints. 

There are remaining issues in developing theoretical frameworks for probing 
gravitational physics with upcoming cosmological surveys and exploiting the constraints 
expected to be obtained. The questions that we need to address include: 

• What theoretical models should we be looking at? 

• Is there a better, more physically meaningful parameterisation of modihcations to 

GR on cosmological scales? 

• How do we include non-linear physics either numerically or analytically? 

• How do we combine various cosmological observations to test gravity? 

• How limited are we to using linear scales? 

• Can we develop novel tests of gravity on non-linear scales? 

As we have seen in this review, we have started to answer these questions and 
have already obtained a number of new constraints on gravity on astrophysical and 
cosmological scales in the last ten years. The next ten years will be an exciting time 
for cosmological tests of gravity thanks to the on-going and future surveys. We may be 
able to detect the breakdown of GR that is responsible for the later acceleration of the 
Universe or GR will be conhrmed as the theory of gravity on cosmological scales. 
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